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ABSTRACT 


The  Campbell  Flats  anticlinorium,  in  the  hanging  wall 
of  tne  Mason  thrust,  underlies  Mount  Hamell  and  Grande 
Mountain,  and  crosses  the  Smoky  River  some  1 5  km  north- 
northwest  of  the  town  of  Grande  Cache,  An  area  of  about  60 
sg  km  was  mapped  in  detail,  field  data  being  stored, 
retrieved  and  processed  using  the  computer.  Existing 
computer-based  procedures  were  used  for  example  to  plot 
maps,  prepare  orientation  diagrams,  divide  the  area  into  22 
domains  within  which  folding  can  be  considered  cylindrical, 
calculate  fold-axes,  rotate  domains,  prepare  cross-sections 
by  downplunge  projection,  and  calculate  stratigraphic 
thicknesses. 

A  computer-based  method  was  developed  to  prepare 
structure-contour  maps  of  coal  seams  or  other  stratigraphic 
horizons  using  only  orientational,  stratigraphic  and 
positional  outcrop  data. 

It  has  four  steps:  (a)  establishment  of  domains  within 
which  the  horizon  is  approximately  cylindrically  folded;  (b) 
construction  of  the  horizon  profile  in  each  domain;  (c) 
projection  of  each  profile  parallel  to  its  fold-axis  in 
order  to  generate  the  coordinates  of  a  set  of  points  on  the 
horizon  in  each  domain;  (d)  computer-contouring  of  the 
resulting  elevations.  The  first  and  fourth  steps  use  known 
procedures.  The  second  step  can  be  carried  out  graphically 
using  a  computer  plot  that  shows  the  geographic  location. 
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stratigraphic  position  and  bedding  trace  of  each  outcrop 
projected  onto  a  plane  normal  to  the  fold-axis. 

Alternatively ,  where  the  stratigraphic  positions  of  the 
outcrops  are  known  precisely  enough  and  the  structure  is  not 
too  complex,  the  computer  can  be  instructed  to  interpolate 
the  horizon  between  the  various  projected  outcrops.  The 
third  step  involves  using  the  appropriate  digitized  profile 
and  fold-axis  to  predict  the  horizon’s  depth  beneath  each 
outcrop  and  saving  these  values  for  the  contouring  stage.  In 
the  case  of  coal  seams,  overburden  thickness  maps  can  also 
be  prepared  using  this  method. 

The  Campbell  Flats  anticlinorium  exposes  strata  ranging 
from  the  Upper  Jurassic  -  Lower  Cretaceous  Nikanassin 
Formation  to  the  Upper  Cretaceous  Kaskapau  Formation.  The 
coal-bearing  Luscar  Formation  was  examined  in  detail  and 
found  to  consist  of  four  units:  (A)  90  m  of  interbedded 
sandstone,  siltstone  and  shale  with  several  thin  coal  seams, 
the  most  important  being  1.5  m  thick;  (B)  30  m  of  marine 
shale,  equivalent  to  the  Moosebar  Formation  of  northeastern 
British  Columbia;  (C)  270  m  of  interbedded  sandstone, 
siltstone  and  shale  with  three  mineable  seams  totalling  8  m 
in  thickness;  (D)  170  m  of  sandstone  with  some  interbedded 

shale  and  conglomerate. 

From  the  structural  standpoint,  folds  within  the 
Campbell  Flats  anticlinorium  are  characterized  by  planar 
limbs,  narrow  hinge  zones  and  some  intralimb  faults.  The 
orientations  of  the  fold-axes  and  axial  planes  of  the  folds 
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change  rapidly  from  one  area  to  another.  These  geometric 
characteristics  and  the  association  with  considerable 


bedding-plane  slip  suggest  that  the  dominan 
mechanism  was  kinking  rather  than  buckling. 
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INTRODUCTION 


Abundant  reserves  of  high  quality  metallurgical  coal 
are  present  in  the  lower  Cretaceous  Luscar  Formation  of  the 
Rocky  Mountain  Foothills  in  west-central  Alberta.  This  coal 
has  been  mined  since  the  beginning  of  the  century  and  will 
undoubtedly  be  exploited  for  many  years  to  come.  Exploration 
and  mapping  for  coal  in  the  Smoky  River  region  began  early 
in  the  1900's  but  because  of  poor  markets  and  the  remoteness 
of  the  area  activity  was  suspended  until  recent  times. 

This  thesis  describes  the  detailed  structure  and 
stratigraphy  of  the  Luscar  Formation,  about  which  little  is 
known,  in  a  small  area  straddling  the  Smoky  River  near 
Grande  Cache  (Fig.  1) .  The  aims  of  the  thesis  were 
essentially  four-fold:  (1)  to  determine  the  stratigraphy  of 
the  Luscar  Formation,  in  particular  the  number  and 
distribution  of  coal  seams;  (2)  to  describe  and  analyse  the 
structure  of  the  Luscar  and  adjacent  formations;  (3)  to 
develop  a  numerical,  computer-based  method  for  preparing 
structure  contour  and  overburden  thickness  maps;  and  (4)  to 
demonstrate  the  usefulness  of  this  and  other  numerical, 
computer-based  methods  in  determining,  displaying  and 
analysing  the  detailed  structure  of  coal  measures.  Field 
work  for  this  study  was  carried  out  in  1977  during  the  four 
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summer  months,  after  which  ten  months  were  spent  in  Edmonton 
as  a  full-time  graduate  student  at  the  University  of 
Alberta, 

The  area  studied  is  near  the  town  of  Grande  Cache, 
about  400  km  w est-northwest  of  Edmonton.  It  is  readily 
accessible  ry  Highway  40  from  Hinton  and  there  is  a  gravel 
airstrip  at  Grande  Cache  suitable  for  light  planes.  The 
Alberta  Resources  Railway  from  Hinton  to  Grande  Prairie 
crosses  the  northwestern  part  of  the  area.  Access  to  the 
northern  part  of  tne  area  is  provided  by  the  road  from 
Grande  Cache  to  the  McIntyre  Mines  plant  on  the  Smoky  River, 
whereas  Highway  40  provides  access  to  the  southern  part. 
Within  the  60  sg  km  of  the  area  itself  there  are  several 
exploration  roads  suitable  only  for  trail  bikes  or  four- 
wheel  drive  vehicles. 

Dominating  the  area  are  Mount  Hamell  (2129  m)  just  west 
of  the  study  area  and  Grande  Mountain  (1987  m)  in  the  south. 
Topographic  relief  is  greatest  and  bedrock  exposure  best  on 
the  slopes  of  the  valley  of  the  Smoky  River  which  flows 
northeastwards  between  the  two  mountains,  and  whose 
elevation  is  as  low  as  930  m.  The  topography  of  the  sides  of 
the  Smoky  River  valley  is  rugged  with  narrow  tributary 
valleys  and  steep-sloped  intervening  spurs.  Beyond  the 
watersheds  enclosing  the  Smoky  River  the  relief  is  generally 
more  subdued,  the  vegetation  much  denser,  and  outcrops  more 
scattered.  Vegetation  is  also  controlled  by  elevation.  Below 
1000  m,  trembling  aspen  are  indicative  of  an  immature 
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montane  rather  than  a  coniferous  montane  forest  which  is 
normally  characteristic  of  these  elevations.  Small  areas  of 
muskeg,  vegetated  by  black  spruce,  willow,  birch  and  an 
assortment  of  grasses,  are  present  where  drainage  is  poor. 

As  the  elevation  increases  to  1500  m  the  vegetation  consists 
increasingly  of  species  typical  to  the  montane  forest 
ecological  zone,  namely  an  assortment  of  conifers,  shrubs, 
herbs  and  grasses.  Above  the  montane  forest,  sub-alpine  to 
alpine  vegetation  is  predominant.  This  zone  is  characterized 
by  small  conifers  and  by  some  shrubs  and  grasses.  Alpine 
meadows  are  common  on  the  flatter  ridge  tops.  The  conifers 
end  at  appoximately  1900  m  and  above  this  the  sensitive 
alpine  tundra  environment  is  encountered  intermixed  with 
areas  of  barren  rock. 

The  region  has  already  been  mapped  on  a  reconnaissance 
scale  by  the  Geological  Survey  of  Canada  (MacVicar, 
1917,1920,1924;  Irish,  1950,1952,1957,1965;  Thorsteinsson , 
1952;  Stott,  1960,1968,1973).  MacVicar  (1917,1920,1924)  and 
Irish  (1965)  provided  brief  accounts  of  the  history  of  the 
area.  Other  studies  that  include  the  Grande  Cache  area  were 
published  by  McEvoy  (1925)  and  McLean  (1977).  Coal  leases 
for  the  area  are  held  by  the  Coal  Division  of  McIntyre  Mines 
Ltd.  who  conducted  exploration  work  during  the  late  1950's 
and  early  1960*s,  and  also  from  1969  to  1975  (Landes,  1962; 
Brown,  1971).  A  number  of  boreholes  and  adits  were  completed 
as  part  of  this  exploration. 


. 
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D AIA  COLLECTION ,  RETRIEVAL  AND  PROCESSING 


Field  Data 


Structural  and  stratigraphic  data  were  collected  in  the 
field  by  finding  outcrops,  plotting  their  positions  on 
aerial  photographs  and  recording  various  observations  made 
at  each  outcrop.  In  addition  certain  stratigraphic  horizons 
were  followed  on  the  ground  and  with  the  aid  of  stereoscopic 
pairs  of  photographs  and  their  traces  marked  on  enlarged 
aerial  photographs.  The  surface  traces  of  formation 
boundaries,  recognizable  stratigraphic  horizons,  faults  and 
axial  planes  of  folds  were  mapped  on  the  photos  and  then 
transferred  to  the  base  map.  These  features  were  continually 
modified  and  updated  as  more  information  became  available 
and  knowledge  of  the  structure  and  stratigraphy  of  the  area 
improved. 

Access  to  tne  outcrops  was  largely  by  foot  along 
traverses  beginning  on  the  main  road  through  the  Smoky  River 
valley.  Midway  through  the  field  season  a  Honda  90  trail 
bike  and  a  Ford  4-wheel  drive  truck  enabled  the  more  remote 
areas  to  be  reached  using  exploration  roads. 
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Aerial  photographs  were  obtained  from  the  Alberta 
Department  of  Energy  and  Natural  Resources.  The  original 
photographs,  taken  in  1970  and  printed  at  a  scale  of 
1:18000,  provided  stereoscopic  coverage,  whereas 
enlargements  at  a  scale  of  1:9000  were  used  in  the  field. 
Topographic  maps  on  a  scale  of  1:6000  with  a  contour 
interval  of  10  ft  and  a  1000  ft  rectilinear  north-south, 
east-west  reference  grid  were  prepared  by  Western 
Photogrammetry  Ltd.  in  Edmonton. 

Throughout  this  study  considerable  emphasis  was  placed 
on  using  compu ter-based  techniques  to  store,  retrieve  and 
analyse  structural  and  other  data.  From  the  outset  field 
data  were  collected  following  a  systematic  and  concise 
procedure  that  was  applicable  to  every  outcrop  visited.  To 
facilitate  collection  and  computerization,  as  many  field 
data  as  possible  were  recorded  in  numeric  code  on  80-space 
sheets  {Fig.  2) . 

Outcrops  were  numbered  consecutively,  each  number  being 
entered  in  columns  1-4.  Altogether  800  outcrops  were 
visited.  Columns  5-8,  designed  for  borehole  identification, 
were  not  used  during  field  data  collection.  A  binary  code 
indicating  (1)  the  formation  and  if  possible  (2)  the 
recognizable  stratigraphic  interval  within  the  formation  was 
entered  in  columns  9-10  (Appendix  1).  Where  possible  the 
stratigraphic  distance  in  feet  of  the  exposed  horizon  above 
the  base  of  the  formation  to  which  it  belongs  was  entered  in 
columns  11-14.  Each  outcrop  was  located  as  precisely  as 
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BEDDING  ORIENTATIONS 
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FIGURE  2.  Field  data  sheet  for  collecting  outcrop 
information. 
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possible  and  then  marked  by  a  pinhole  and  the  location 
identified  on  the  back  of  the  appropriate  1:9000  photo. 

After  each  outcrop  location  had  been  transferred  to  the 
topographic  base  map,  its  coordinates  in  feet  were  obtained 
and  recorded  in  columns  15-29  of  the  appropriate  field  data 
sheet.  The  X  and  Y  coordinates,  referred  to  horizontal  axes 
pointing  due  east  and  due  north,  respectively,  with  an 
origin  at  the  UTH  grid  point  19,600,000N  and  150,000E  of  the 
120th  meridian,  were  found  with  the  aid  of  the  rectilinear 
grid,  a  Douglas  protractor  and  an  engineer’s  scale.  The  Z 
coordinate,  in  feet  above  sea-level,  was  measured  by 
interpolating  between  contours.  Errors  in  the  coordinates  of 
an  outcrop  stemmed  from  four  possible  sources: 

1.  Locating  the  outcrop  incorrectly  on  the  aerial  photo. 

2.  Transferring  the  outcrop  incorrectly  from  the  aerial 
photo  to  the  topographic  map. 

3.  Reading  the  coordinates  incorrectly  from  the  topo¬ 
graphic  map. 

4.  Inaccuracies  in  the  topographic  map. 

The  first  three  types  of  error  combined  normally  should  not 
be  more  than  33  ft  (10  m) ;  the  fourth  is  unlikely  to  exceed 
6  ft  (2  m) .  The  likely  error  in  the  position  of  an  outcrop 
as  given  by  its  coordinates  was  estimated  in  feet,  divided 
by  ten  and  placed  in  column  30. 

At  each  outcrop  up  to  six  bedding  plane  orientations 
were  measured  within  an  area  of  1 0  sg  m  and  recorded  in 
columns  15-29  of  the  data  sheet.  Planar  orientation  data 
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were  measured  as  numerical  dip-direction  and  dip  with  a 
Ereinerger  compass  scaled  in  grads  rather  than  degrees. 
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The  presence  of  additional  information  recorded  in  a 
conventional  field  notebook  was  indicated  on  the  field  data 
sheet  by  a  1=yes  or  G=no  entry  in  columns  31-38  (Appendix 
1) .  Such  information  concerned  (a)  the  orientations  of 
bedding  plane  linears,  mesoscopic  folds,  faults  and  joints, 
(b)  lithologic  descriptions,  (c)  fossils,  (d)  general 
observations  and  (e)  photographs.  Where  bedding  plane 
linears  were  observed,  the  orientations  were  recorded  as 
multiple  readings  of  the  pitches  of  these  linears  measured 
clockwise  on  the  bedding  plane.  At  outcrops  displaying 
mesoscopic  folds,  as  wide  a  range  as  possible  of  bedding 
orientations  was  recorded  along  with  the  orientation  of  the 
axial  plane.  Multiple  readings  of  fault  and  joint  plane 
orientations  were  recorded  as  numerical  dip-direction  and 
dip  along  with  the  pitches  of  any  associated  slickensides. 
Repetitive  cleat  plane  orientations  were  obtained  at  some  of 
the  coal  outcrops.  Both  35  mm  slides  and  Kodak  EK-6  instant 
pictures  were  taken  in  the  field.  The  EK-6  pictures  were 
extremely  useful  for  drawing  on  and  noting  field 
observations,  but  the  camera  itself  was  cumbersome  and  the 
pictures  were  generally  of  poor  quality. 

The  data  sheets  were  very  convenient  to  use  but  could 
be  improved  by'  a  few  minor  modifications.  The  borehole 
identification  could  be  incorporated  into  the  station  number 
entry  and  the  spaces  between  the  bedding  orientations 
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eliminated-  In  place  of  these  the  date  and  codes  indicating 
the  weather  conditions,  time  of  day  and  mesoscopic  fold 
geometry  could  be  included-  Fold  information  could  be  sub¬ 
divided  into  four  separate  entries:  (a)  type  of  folding  - 
anticline,  syncline,  monocline  or  series  of  folds;  (b) 
orientation  of  the  fold  -  upright  horizontal,  upright 
plunging,  vertical,  inclined  horizontal,  inclined  plunging, 
recumbent  or  reclined;  (c)  style  of  folded  surfaces  - 
chevron,  box  or  rounded;  and,  (d)  style  of  folded  layers  - 
parallel  or  similar-  This  would  provide  the  geologist  with  a 
very  accessible  index  on  the  geometry  of  mesoscopic  folds. 

It  would  also  be  possible  to  examine  the  productivity  of  the 
geologist  under  various  weather  conditions  and  times  of  the 
day  and  monitor  changes  in  productivity  throughout  the  field 
season. 


Borehole  Data 


Borehole  data  in  the  thesis  area  are  few  and  generally 
of  poor  quality-  Although  6  of  the  49  boreholes  examined  had 
gamma  ray  and  density  logs,  only  elevations  and  written 
descriptions  of  the  intersected  coal  seams  were  available  in 
the  case  of  the  remaining  43  boreholes.  Borehole  deviation 
was  generally  taken  into  account  in  calculating  the  X,  Y  and 
Z  coordinates  of  the  coal  intersections.  However,  10  of 
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these  boreholes  drilled  on  Grande  Mountain  were  not  provided 
with  deviation  records.  Although  these  holes  were  assumed  to 
be  vertical,  with  bedding  dipping  approximately  20  NE  and 
the  drillholes  being  up  to  300  m  long,  some  deviation  from 
this  orientation  is  to  be  expected. 

During  previous  exploration,  approximately  10  adits 
were  driven  for  bulk  sampling  of  the  coal.  Eor  this  study 
the  adits  served  only  as  artificial  outcrops  of  the  coal 
seams. 


Construction  and  Preliminary  Processing  of  Data  Files 


The  first  step  upon  arrival  back  at  the  university  was 
to  construct  a  master  data  file  containing  all  the  codeable 
and  numeric  information  in  an  accessible  form  for  retrieval 
and  processing.  A  master  file  DATA  was  constructed  by 
reading  computer  punch  cards  prepared  from  the  800  outcrop 
data  sheets.  The  format  of  all  but  the  first  and  last  lines 
of  DATA  is  identical  to  the  field  data  sheets  with  column 
numbers  in  the  file  equivalent  to  spaces  on  the  data  sheets 
(see  Appendix  1). 

The  storage  file  BOREHOLE,  containing  important 
stratigraphic  intersections  and  the  borehole  surface 
locations  was  constructed  from  the  borehole  information  (see 
Appendix  2) .  The  coordinates  of  the  intersections  of  faults 
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and  the  tops  of  coal  seams  and  the  Caaomin  Formation  were 
calculated  from  oorehole  surface  coordinates,  borehole 
orientation  and  deviation  and  drilled  depths. 

The  field  data  contained  in  the  file  DATA  were 
partially  processed  and  stored  in  the  file  MEANN  whose 
format  was  acceptable  as  input  to  the  domain  analysis 
computer  packages  available  in  the  Department  of  Geology. 

The  file  DATA  was  first  sorted,  by  the  systems  package 
*SORT*,  according  to  stratigraphic  code  and  was  input  to  the 
program  DEIET.  Using  a  series  of  temporary  files  and  the 
programs  SORT,  MEANS  and  RESORT,  the  file  MEANN  was  then 
constructed  (Fig.  3). 

The  program  DELET  (Appendix  3)  converts  the  orientation 
data  measured  in  grads  to  degrees,  assigns  a  letter  to  each 
stratigr aphic  interval  and  places  the  output  in  a  temporary 
file  -DATA1.  The  program  SORT  (Appendix  3)  using  -DATA  1  as 
input  constructs  a  temporary  file  -DATA2  containing  the 
repetitive  bedding  orientations  at  each  outcrop  in  a  format 
acceptable  as  input  to  the  program  MEANS  (available  in  the 
Department  of  Geology) .  A  mean  bedding  orientation  for  each 
outcrop  is  calculated  by  this  program  through  a  process  of 
vector  addition  of  the  direction  cosines  of  the  poles  to 
bedding  and  correcting  the  resultant  mean  vector  to  unity. 
Fisher's  concentration  parameter  K,  which  measures  the 
scatter  of  the  bedding  orientations  measured  at  each 
outcrop,  is  also  calculated  (Mardia,  1972,  p.  251).  K  values 
calculated  from  outcrop  data  in  the  thesis  area  ranged  from 
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FIGURE  3  I  FLGH  OI/tOR.W  OF  FViOCEOlMtS  FOR  DATA  PROCESSING 
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about  30  (indicating  a  scatter  in  orientations)  to  over 
40,000  (suggesting  very  little  difference  in  individual 
bedding  orientations) .  The  output  from  the  program  MEANS  (in 
the  temporary  files  -MEAN  and  -KV)  along  with  certain  data 
from  the  file  -DATA1  are  then  rearranged  by  the  program 
RESORT  (Appendix  3)  and  the  output  placed  in  the  file  MEANN 
whose  format  is  given  in  Appendix  4. 

Geological  maps  (e.g.  Fig.  4,  in  pocket)  showing  the 
locations  of  outcrops,  mean  bedding  orientations  and,  as  an 
option,  outcrop  numbers,  were  produced  on  the  Calcomp 
Plotter  or  Tectronix  Graphics  Terminal  using  the  program 
DRAW  (see  Appendix  5)  ,  Input  data  for  this  program  were 
obtained  from  the  file  MEANN.  Similar  maps  showing  borehole 
locations,  bedding  plane  striae  and  mesoscopic  fold-axes 
were  also  produced  by  DRAW.  The  format  of  the  input  data  is 
the  same  in  each  case  except  that  the  dip-directions  and 
dips  are  replaced  by  zeros  for  boreholes  or  by  trends  and 
plunges  for  striae  and  fold-axes. 


Data  Retrieval 


Two  methods  were  used  to  retrieve  data  from  the  file 
MEANN  and  rearrange  them  in  a  form  amenable  tc  further 
processing  and  analysis.  The  first  method  involved  using  the 
program  SELECT  (see  Appendix  3)  to  retrieve  information  by 
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outcrop  number.  One  of  the  input  files  contained  the  desired 
outcrop  numbers  and  the  second  the  data  from  which 
information  was  to  be  retrieved.  The  format  of  SELECT  had  to 
be  modified  as  different  data  were  required  from  different 
files.  The  second  method  of  data  retrieval  was  by  means  of 
the  file  editor  and  copy  commands.  In  most  instances  the 
data  files  were  sorted  according  to  some  parameter  before 
editing  or  copying  took  place.  During  later  stages  of 
structural  analysis,  the  domain  files  were  edited  and  the 
appropriate  borehole  information  was  inserted. 


Data  Processing 


Determination  of  Cylindricity 
Using  the  geological  map  and  a  qualitative  knowledge  of 
the  structural  geology  of  the  area,  domains  were  delineated 
within  which  folding  was  expected  to  be  cylindrical. 

Separate  domain  files  whose  format  was  identical  to  MEANN 
and  BOREHOLE  and  containing  all  the  outcrop  and  borehole 
data  for  each  domain  were  constructed  using  both  data 
retrieval  techniques  and  then  analysed  for  cylindrical 
folding. 

Individual  domains  were  processed  initially  by  the 
program  CYLTEST,  available  in  the  Department  of  Geology, 
which  determines  a  best-fit  fold-axis  and  some  statistical 
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parameters  that  are  used  in  testing  for  cylindricity.  The 

vector  of  direction  cosines  of  the  pole  to  bedding  is 

determined  for  each  outcrop  in  a  domain.  The  column  vector 

of  direction  cosines  for  each  bedding  pole  is  then 

postmultiplied  by  the  row  vector  and  the  products  are  summed 

» 

to  form  a  symmetric  3x3  matrix  which  has  three  eigen  values 

and  associated  eigen  vectors.  If  <£.•  is  the  angle  between  the 

ith  of  p  poles  to  bedding  and  an  estimated  fold-axis,  the 

n 

best  fit  fold-axis  is  obtained  when  £  coszcf>.is  minimized. 

1  =  1  1 

This  minimized  sum  of  squares  of  cos2^-  is  equal  to  the 
minimum  eigen  value  { A3  )  and  the  associated  3^element  eigen 
vector  contains  the  direction  cosines  of  the  best-fit  fold- 
axis  (e.g.  Cruder,  1968,  p.  143-148). 

For  a  perfectly  cylindrical  fold,  with  no  scatter  of 
bedding  poles  about  the  plane  normal  to  the  fold-axis,  the 
value  of  A3  is  0.  Because  of  measurement  error  and  bedding 
plane  roughness  on  a  scale  less  than  that  of  the  folding 
being  investigated,  this  ideal  situation  is  unlikely  to  be 
encountered,  even  in  the  case  of  cylindrical  folds.  Two 
statistical  tests  are  used  to  determine  the  likelihood  of 
the  observed  scatter  of  bedding  poles  being  due  to  non- 
cylindricity  (e.g.  Charlesworth  et  al,  1976). 

1.  The  coplanarity  test  determines  the  likelihood  of  the 
deviation  from  cylindricity  being  due  to  measurement 
error  and  bedding  plane  roughness  on  a  scale  less  than 
that  of  an  individual  outcrop.  Fisher's  concentration 
parameter  (K)  and  derived  from  single  randomly 
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selected  original  bedding  orientations,  1  per  outcrop, 
are  required  to  determine  a  chi-square  statistic  used  to 
test  the  null  hypothesis  of  cylindrical  folding. 

2.  The  coaxiality  test  determines  the  likelihood  of  the 
deviation  from  cylindricity  being  due  to  measurement 
error  and  bedding  plane  roughness  on  a  scale  larger  than 
that  of  an  outcrop.  Each  domain  is  divided  into  two 
segments  (a  and  b)  with  p^  and  outcrops  and  two 
values  of  A3  (  A3a  and  A3^)  and  the  respective  fold-axes 
are  determined.  To  accept  the  cylindrical  model  two 
conditions  must  be  met:  first,  the  scatter  of  poles  in 
each  segment  about  planes  normal  to  the  respective  fold- 
axes  must  be  equal  and,  second,  the  deviation  from 
cylindricity  for  the  entire  domain  must  be  accounted  for 
by  the  deviation  from  cylindricity  in  the  two  segments, 
for  practical  purposes  the  results  of  these  two  tests 
are  useful  but  should  not  be  rigorously  adhered  to  when 
accepting  or  rejecting  the  cylindrical  model.  Other  criteria 
include  (a)  the  relative  values  of  Ai  /  A2  and  \3 /  (b)  the 
value  of  the  standard  scattering  angle  -  the  standard 
deviation  of  the  angles  between  the  bedding  poles  and  the 
plane  normal  to  the  fold-axis,  and  (c)  the  value  of  the 
angle  between  the  fold-axes  of  different  segments  of  the 
domain  (Kilby,  1978,  p.  26-29)  .  Throughout  this  study 
cylindricity  was  assumed  when: 

1.  A 2  was  less  than  an  order  of  magnitude  smaller  than  \l 
and  more  than  an  order  of  magnitude  larger  than  A3  ; 
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2.  the  standard  scattering  angle  was  less  than  6°; 

3.  the  angle  between  the  fold-axes  of  separate  segments 
of  a  domain  was  less  than  5°;  and 

4.  the  null  hypothesis  of  cylindricity  associated  with 
one  of  the  statistical  tests  could  not  be  rejected. 


Profiles  and  Composite  Sections 
The  program  SECTS,  available  in  the  Department  of 
Geology,  was  used  to  obtain  plots  showing  for  each  domain 
down  plunge  projections  of  outcrops  and  borehole 
intersections  within  the  domain  (see  Appendix  6)  .  Input  to 
SECTS  consists  of  the  domain  file  containing  station 
locations  and  bedding  orient atioD s,  the  fold-axis 
orientation,  the  normal  to  the  plane  of  projection  and  a 
scaling  factor  to  determine  the  size  of  the  plot.  These 
plots  can  then  be  used  to  construct  profiles  and  sections  of 
any  orientation  (Fig.  5) .  For  each  data  observation  point  in 
the  input  file,  the  program  SECTS  transforms  the  map 
coordinates  to  profile  coordinates  (X  and  Y)  and  calculates 
the  pitcn  of  bedding  on  the  profile.  At  each  location  on  the 
plane  of  the  profile  a  short  line  to  represent  the  pitch  of 
bedding  and  an  identification  symbol,  present  in  column  1  of 
the  input  file,  are  drawn  by  the  Calcomp  Plotter  or 
Tectronix  Graphics  Terminal.  Locations  without  bedding 
orientations  are  represented  by  a  cross  rather  than  a  line. 
Additional  X,  Y  and  Z  coordinates  of  coal  seams  were  read 
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□  OMR  I N  22 


FIGURE  5.  Downplunge  projection  of  domain  22  with  Luscar 
Units  A,  B,  and  C,  and  #4  and  #10  coal  seams. 
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from  the  geological  map  and  the  locations  plotted  as 
outcrops  without  orientations. 

Composite  plots  for  several  domains  were  ultimately 
constructed  using  the  programs  DOMROT,  available  in  the 
Department  of  Geology,  and  SECT5.  The  domains  to  be 
represented  were  first  rotated  using  DOMROT  so  that  their 
fold-axes  became  parallel.  The  various  rotated  outcrops  and 
borehole  inter sections  were  then  projected  onto  a  single 
plane  using  SECT5.  Input  to  DOMROT  consists  of  the  domain 
file  containing  the  data  to  be  rotated,  the  present  and 
•’new*'  fold-axes  and  the  coordinates  of  the  point  about  which 
the  rotation  is  to  be  performed.  The  output  file  contains 
the  rotated  coordinates  and  bedding  orientations. 

Pi  Diagrams 

The  standard  files  for  domain  analysis  are  input  to  the 
program  STEREO  which  returns  the  temporary  file  -DATA 
consisting  of  a  title,  an  identification  of  data  type  and 
the  orientation  data.  The  program  FINPUT  using  -DATA  as 
input  returns  the  direction  cosines  of  the  orientation  data. 
The  program  PLOT  produces  a  grid  of  numbers,  each  of  which 
is  the  equal-area  projection  of  the  density  of  bedding  poles 
at  a  certain  counting  location  on  the  surface  of  the 
reference  hemisphere.  The  diagram  is  completed  by  contouring 
the  grid  and  enclosing  it  in  a  circle  whose  diameter  is  8  in 
(see  Kilby,  1978,  p.  34). 
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Fold  Geometry 

The  fold-axis  orientation  in  each  domain  can  be 
determined  in  one  of  two  ways:  (a)  by  using  the  program 
CYLTEST,  or  (b)  by  using  the  program  STEEEO  and  the  source 
programs  f INPUT  and  EIGEN,  both  of  which  use  an  eigen  value 
routine  (p.  16) .  The  best-fit  fold-axis  and  the  pitch  of  the 
axial  plane  measured  on  the  profile  normal  to  the  fold-axis 
(p.  18),  are  used  to  define  the  axial  plane  orientation.  The 
apical  angle  can  be  measured  and  the  limb  shape  determined 
from  the  down  plunge  projection. 

Stratigraphic  Thicknesses 

The  poor  exposure  in  parts  of  the  thesis  area  made  the 
stratigraphic  position  of  certain  outcrops  difficult  to 
determine.  The  API  program  THICK,  implementing  an  automated 
procedure  of  the  Busk  technique,  was  used  to  calculate 
approximate  stratigraphic  thicknesses  between  pairs  of 
exposures  (Kilby,  1978,  p.  34^35) .  The  difference  between 
the  radii  of  cylinders  with  a  common  axis  at  the 
intersection  of  planes  normal  to  the  two  bedding  planes,  and 
with  surfaces  passing  through  the  outcrops,  is  the  interval 
thickness.  In  situations  where  the  angle  between  the  two 
bedding  planes  is  less  than  5°,  the  thickness  is  taken  as 
the  perpendicular  distance  between  the  planes  passing 
through  the  two  outcrops  whose  orientation  is  the  mean  of 
that  at  these  outcrops.  This  method  assumes  that  folding  is 
concentric  so  in  this  map  area,  where  the  predominant  styles 
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of  folding  are  chevron  and  rounded  chevron  with  narrow  hinge 
zones,  both  exposures  must  be  located  on  the  same  limb  for 
this  technique  to  be  applicable. 

Mesoscopic  Structural  Processing 
Mesoscopic  information  concerning  bedding  plane  striae, 
mesoscopic  folds,  mesoscopic  faults,  joints  and  the  Cadomin- 
Nikanassin  unconformity  were  collected  and  subsequently 
processed  and  analysed.  Approximately  300  orientations  of 
striae  on  bedding  and  fault  planes  were  collected  at  25 
outcrops.  Temporary  files  containing  a  title  in  line  1,  the 
bedding  or  fault  plane  orientation  in  the  first  6  columns  of 
line  2  and  the  pitches  of  the  striae  in  the  first  3  columns 
of  the  remaining  lines  were  constructed.  Each  of  these 
temporary  files  were  processed  by  the  program  SLICKS  (see 
Appendix  3)  with  the  output  going  to  the  file  SLICKENS.  The 
output  consists  of  a  title  in  line  1,  the  data  type  in  line 
2  (in  this  case  IYPE=1)  and  the  trend  and  plunge  of  the 
striae  in  the  first  6  columns  of  all  the  remaining  lines 
except  the  last  in  which  an  end  of  data  statement  is 
present.  The  file  SLICKENS  was  then  copied  to  -DATA  and 
processed  by  the  source  programs  EINPUT  and  EIGEN.  The  file 
MESODATA  was  constructed  containing  the  X,  Y  and  Z 
coordinates  in  columns  2-16,  the  mean  trend  and  plunge  of 
the  linears  at  each  outcrop  in  columns  18-23  and  a  2  in 
column  80.  The  trends  and  plunges  of  the  striae  were  plotted 
on  a  mesoscopic  data  map  by  the  program  DRAW  (fig.  6). 
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FIGURE  6 
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Mesoscopic  fold  data  were  of  two  types:  (a)  mesoscopic 
fold-axis  orientations,  and  (b)  bedding  orientations 
measured  across  the  folds-  Multiple  fold-axis  orientations 
were  placed  in  temporary  files  and  processed  as  linears  by 
the  source  programs  FINPUT  and  EIGEN  to  determine  a  mean 
fold-axis  orientation  for  each  outcrop  with  mesoscopic 
folds-  Each  set  of  bedding  orientations  was  placed  in  a 
temporary  file  and  processed  by  the  programs  FINPUT  and 
EIGEN  to  determine  a  best  fit  fold-axis  orientation.  The 
fold-axis  orientations  were  put  into  the  file  MESODATA  with 
the  same  format  as  the  striae  data  but  with  a  3  in  column 
80. 

Multiple  fault  planes  were  processed  by  FINPUT  and 
EIGEN  to  determine  a  mean  orientation  at  each  outcrop  where 
faulting  was  observed.  The  angularity  between  bedding  and 
the  fault  plane  was  also  determined  at  these  outcrops. 

A  mean  joint  plane  orientation  was  obtained  from  FINPUT 
and  EIGEN  for  each  joint  set  observed  and  a  pi  diagram  of 
the  joint  plane  orientations  was  produced  using  the  program 
PLOT. 

Computer  programs  written  specifically  for  analyses  of 
structural  data  included  in  this  thesis  are  listed  with 
appropriate  run  commands  in  Appendix  4. 
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SIR  UCTURE  CONTOURING 


Contour  maps  showing  the  configuration  of  a  coal  seam 
and  the  variation  of  its  depth  below  the  topographic  surface 
are  used  extensively  in  coal  exploration.  A  major  aim  of 
this  thesis  was  to  develop  a  numerical,  computer-based 
technique  for  constructing  such  maps  from  an  input  of 
orientation,  positional  and  stratigraphic  data.  Most  contour 
maps  are  constructed,  either  manually  or  numerically  (e.g. 
Gold,  1978) ,  from  irregular  sets  of  points  in  two- 
dimensional  space,  at  each  of  which  the  numerical  value  of 
the  quantity  to  be  contoured  is  known.  However,  the  above 
contouring  method  can  be  implemented  in  the  case  of  coal 
seams  only  after  extensive  drilling.  What  is  needed  is  a 
method  that  uses  as  input  an  irregular  set  of  outcrops  and 
borehole  inter sections  of  known  geographic  location, 
stratigraphic  position  and,  in  the  case  of  outcrops,  bedding 
orientation . 

In  the  computer-based  method  of  structure  contouring 
from  outcrop  and  borehole  data  developed  by  Robinson  (1976), 
coordinates  on  the  horizon  to  be  contoured  are  estimated  by 
simulation  of  the  known  fold  geometry  within  a  cylindrical 
domain.  The  method  developed  in  this  thesis  has  four  steps: 
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1.  division  of  the  area  into  domains  within  which  the 
horizon  to  be  contoured  is  cylindrically  folded,  and 
calculation  of  the  fold-axis  in  each  domain; 

2.  construction  of  the  profile  of  the  horizon  to  be 
contoured  in  each  domain; 

3.  projection  of  each  profile  parallel  to  its  fold-axis 
in  order  to  generate  the  coordinates  of  a  set  of  points 
on  the  horizon;  and 

4.  computer  contouring  of  the  resulting  elevations. 

The  first  step  has  already  been  dealt  with  (p.  14)  and  the 
fourth  is  already  well  known.  The  second  and  third  steps  are 
described  below. 

The  Profile 


Profiles  may  be  constructed  either  graphically  or  using 
a  numerical,  computer-based  procedure.  The  graphical 
technique  is  used  where  the  exact  stratigraphic  position  of 
most  outcrops  and  borehole  intersections  is  uncertain  or 
where  structural  complexities  require  considerable 
interpretation  on  the  part  of  the  person  drawing  the 
profile.  The  profiles  can  be  drawn  on  computer  plots  that 
show  the  positions  of  each  outcrop  and  borehole  intersection 
projected  parallel  to  the  fold-axis  and  where  possible  the 
traces  of  bedding  (p.  18) .  Using  as  a  guide  the  traces  of 
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bedding  and  the  stratigraphic  positions  of  the  outcrops  and 
borehole  intersections,  a  best-fit  line  representing  the 
trace  of  the  horizon  to  be  contoured  is  drawn  on  each  plot. 
The  resulting  lines  are  then  digitized,  generally  at 
irregular  distances  along  the  horizontal  (X)  axis,  care 
being  taken  to  include  the  hinge  and  inflection  points  of 
any  folds,  and  each  digitized  point  is  given  equal  weight 
(see  below) .  The  pitch  of  the  horizon  at  each  digitized 
point  (X  Y  )  is  estimated  from  the  profile. 

The  numerical  method  for  constructing  profiles  can  be 
used  where  the  exact  stratigraphic  position  of  most  outcrops 
is  known  and  where  the  fold  is  either  similar  or  parallel  in 
style.  It  consists  of  the  following  steps: 

(1)  a  computer  line  file  is  constructed  from  the  profile 

coordinates  of  each  outcrop  and  borehole  intersection  in 
the  domain  projected  normally  onto  the  plane  of  the 
profile  and,  where  possible,  the  pitch  of  the  trace  of 
bedding  and  the  precise  stratigraphic  position.  To  do 
this  each  domain  file  is  processed  by  SECT5  (p.  18)  .  The 

profile  coordinates  and  the  pitches  of  bedding  from  the 
output  file  are  then  combined  with  the  precise 
stratigraphic  positions,  outcrop  numbers  and 
stratigraphic  identifications  from  the  domain  file  by 
the  program  DATAPORM  (see  Appendix  7) . 

(2)  the  profile  coordinates,  bedding  pitch  and  precise 
stratigraphic  position  in  each  line  of  the  above  file 
are  used  to  determine  the  profile  coordinates  (X,.  Y ^  ) 


. 
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of  one  point  on  the  horizon  to  be  contoured.  In  Figure 
7,  E  is  the  normal  projection  of  an  outcrop  onto  the 
plane  of  the  profile,  i.e.  the  XY  plane.  The  coordinates 
of  E  referred  to  the  axes  X  Y  and  the  origin  0  are  (xg 
y  ) .  The  pitch  of  bedding  measured  clockwise  on  the  XY 
plane  is  p  °  and  that  of  the  axial  plane  is  pb°  .  The 
axes  a  and  b  intersecting  at  E  have  pitches  of  pb°  and 
p^+90°,  respectively.  From  Figure  7,  if  the  horizon 
exposed  at  E  is  a  distance  d  stratigraphically  above  the 
horizon  to  be  contoured,  the  coordinates  of  a  point  F  on 
this  horizon,  referred  to  the  axes  a  b  and  to  an  origin 
at  E,  are  (0  d)  in  the  case  of  parallel  folding.  In  the 
case  of  similar  folding  the  point  corresponding  to  F  is 
G,  and  its  coordinates  are  (d-ftan  (p  -p  )  d)  .  Using  well 
known  transformation  equations,  the  coordinates  {x  y  ) 
of  F  and  (x  y  )  of  G  referred  to  the  axes  X  Y  and  to  an 
origin  at  E  are  given  by. 


£  M]  x  [  0  d  ]•  =  [  x  y  ] 1 

and  £M]  [d-Jtan(p  -p  )  d]1  =  [x  y  ]  ■ 


cos{pb)  cos  (180-pb) 

where  [ M ]  = 

cos(270-pb)  cos(pb  +  90) 


The  coordinates  of  F  and  G  referred  to  the  axes  X  Y  and 
to  an  origin  at  0  can  be  obtained  by  adding  x  y  to  x  y 
and  the  pitch  of  bedding  at  F  and  G  is  the  same  as  that 
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FIGURE  7 


Determination  of  the  position  of  the  coal  seam 
from  positional,  orientation  and  stratigraphic  information. 
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at  E.  5nia.ll  changes  to  the  above  equations  are  required 
when  p>  90°.  A  weight  is  calculated  based  on  the 
stratigraphic  interval  d  separating  the  outcrop  or 
borehole  intersection  and  the  horizon  to  be  contoured. 

(3)  the  points  produced  in  step  (2)  display  a  certain 
amount  of  scatter  about  a  best-fit  line  on  the  profile 
representing  the  coal  seam.  Coordinates  and  slopes  of 
points  lying  on  this  best-fit  curve  are  calculated  at 
regular  intervals  along  the  X-axis  by  projecting  the 
profile  coordinates  (X^  Y ^ )  parallel  to  their  slopes  to 
each  X  increment  and  weighting  each  value  relative  to 
the  distance  of  the  projection  along  the  X-axis  and  the 
stratigraphic  weight  calculated  in  step  (2) .  The  slope 
at  a  point  on  the  curve  is  estimated  by  weighting  the 
input  slopes  with  the  stratigraphic  and  projection 
weights . 

The  procedure  of  calculating  the  coordinates  of  a  coal  seam 
on  the  profile  discussed  here  is  by  no  means  a  general  case. 
Anyone  adopting  this  approach  has  to  examine  the  technique 
closely  and,  if  necessary,  modify  it  to  be  consistent  with 
the  geometry  of  the  folds  being  studied.  Plots  shoving  the 
positions  of  points  on  the  coal  seam  are  produced  in  both 
steps  (2)  and  (3) ,  unsmoothed  and  smoothed,  respectively 
(see  program  HORIZON,  Appendix  7)  . 
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The  Three-dimensional  Array 


The  procedure  for  extrapolating  the  profiles  to 
generate  three-dimensional  sets  of  points  representing 
elevations  and  depths  is  outlined  below: 

(1)  a  line  file  containing  the  coordinates  of  a  set  of 
points  (x  y  z  )  on  the  topographic  surface  in  each 
domain  is  built.  This  array  could  form  a  uniform  grid 
and  as  such  would  have  to  be  built  from  scratch  using 
the  topographic  map  of  the  area.  Alternatively  the  array 
could  be  built  using  the  coordinates  of  outcrops  and 
borehole  locations,  supplemented  by  some  points  along 
the  axial  traces  of  folds,  and  as  such  be  irregular.  The 
second  method  was  preferred  not  only  because  it  was  less 
onerous  but  also  because  the  array  of  points,  used 
subsequently  as  input  to  the  contouring  procedure  itself 
(see  below) ,  is  widely  spaced  where  outcrop  and  borehole 
control  is  poor  and  closely  spaced  where  control  is 
good,  as  indeed  should  be  the  case, 
the  next  step  is  tc  calculate  the  elevation  of  the 
horizon  to  be  contoured  directly  beneath  each  of  the 
points  (x  y  z  ) .  To  do  this  the  program  SLOPE  (see 
Appendix  7)  is  used  to  find,  first,  the  profile 
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coordinates  (X  Y  )  of  the  points  (x  y  z  )  projected 
normally  onto  the  plane  of  the  profile,  and  second,  the 
B  (Y)  value  of  the  horizon  to  be  contoured  at  the  node  Q 
.  The  horizon  is  projected  from  each  of  its  profile 
points  (X..  Y  ^ )  (determined  either  graphically  or 
numerically)  parallel  to  the  individual  pitches  to  the 
node  Q  (Fig.  8).  A  weight  is  calculated  for  each 
projection  based  upon  the  distance  of  the  projection 
along  the  X-axis.  A  single  point  (X  B  )  of  the  horizon 
on  the  profile  is  calculated  using  these  projections  and 
the  stratigraphic  and  projection  weights.  This  point  is 
the  normal  projection  of  the  point  [ x  y  z-  (B/cos  p)  ] 
where  p  is  the  plunge  of  the  fold-axis,  so  the  elevation 
and  depth  of  coal  at  the  point  (x  y  z  )  are  (z-E)/cos  p 
and  B/cos  p,  respectively.  The  slope  of  the  horizon 
beneath  the  point  (x  y  z  )  is  determined  by  similarily 
weighting  the  individual  slopes  of  the  profile  points 
(X.j  Y.j)  on  the  coal  seam. 

A  plot  of  the  downplunge  projections  of  the  set  of  points 
located  on  the  surface  of  the  coal  seam  is  produced.  The 
stratigraphic  and  projection  weights  are  calculated  with  the 
equation  Wt=e  (Newton,  1973)  where  A  is  a  constant  and  x 
is  the  distance  between  the  known  and  interpolated 
coordinates.  The  amount  of  smoothing  and  influence  of 
neighbouring  points  can  easily  be  adjusted  by  changing  the 
value  of  A.  The  maximum  possible  weight  using  this  procedure 


is  1.0. 
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Estimation  of  the  profile  coordinates  of  the  coal  seam 
at  any  node  O  from  known  locations  and  orientations. 


FIGURE  8 
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Contouring  the  Three-dimensional  Array 


The  sets  of  coal  seam  elevations  or  depths  calculated 
for  each  domain  were  combined  for  final  contouring  with  a 
triangular  element  routine  provided  by  Dr-  C.  M.  Gold  (Gold 
et  al,  1977) ,  Department  of  Geology,  University  of  Alberta. 
The  procedure  involves  interpolation  of  the  surface  within 
irregular  triangular  domains  constructed  from  the 
coordinates  and  slopes  of  the  horizon  provided  by  the 
program  SLOPE.  Each  individual  array  was  constructed  so  as 
not  to  overlap  but  only  to  border  any  neighbouring  domains 
and  allow  a  smooth  continuum  of  contours  between  cylindrical 
domains  (see  Fig.  14-19). 

These  techniques  can  readily  be  expanded  to  produce 
isopach  and  stripping  ratio  maps.  Slope  calculation  at  each 
map  location  could  be  modified  to  a  weighted  eigen  value 
routine  rather  than  the  weighted  average  presently  being 
used. 
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STRATIGRAPHY 


Some  1250  m  of  marine  and  non-marine  strata  ranging 
from  the  Jurassic  and  Lower  Cretaceous  Nikanassin  formation 
to  the  Upper  Cretaceous  Kaskapau  Formation  are  exposed  in 
the  study  area.  Lithic  sandstone,  siltstone  and  shale  are 
more  common  than  conglomerate  and  coal.  Four  major 
tr ansgressive  and  three  major  regressive  sequences  can  be 
recognized.  The  transgressive  phases  are  represented  by  the 
lower  200  m  of  the  exposed  part  of  the  Nikanassin  Formation, 
a  marine  shale  band  in  the  Luscar  Formation  and  the  marine 
shales  of  the  Shaftesbury  and  Kaskapau  Formations.  The  upper 
160  m  of  the  Nikanassin  Formation,  the  Cadomin  Formation  and 
lower  95  m  of  the  Luscar  Formation  constitute  one 
continental  regime.  The  others  are  represented  by  the  upper 
475  m  of  the  Luscar  Formation  and  by  the  Dunvegan  Formation. 
A  number  of  coal  seams  are  present  within  the  lower  310  m  of 
the  Luscar,  three  of  which  are  of  sufficiently  high  quality 
and  thickness  to  warrant  mining  and  are  currently  worked  by 
McIntyre  Mines  north  of  the  study  area. 
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Nikanassin  Formation 


Some  365  m  of  the  upper  Nikanassin  Formation  are 
exposed  directly  above  the  Grande  Mountain  thrust  and  in  the 
core  of  the  Campbell  Anticlinorium.  Although  no  distinct 
mappable  units  were  encountered,  the  upper  160  m  are  non¬ 
marine  whereas  a  highly  variable  flyschoid  sequence  makes  up 
the  remainder  of  the  formation.  The  main  lithologies 
represented  are  shale,  carbonaceous  shale,  siltstone  and 
lithic  and  argillaceous  sandstone. 

The  non-marine  sequence  contains  carbonaceous  shale  and 
a  few  thin  lenticular  coal  bands.  Its  thin,  cross-bedded 
sandstones  are  commonly  ripple  marked  and  contain  a  few 
thin,  discontinuous  pebbly  bands.  The  marine  sequence  is 
composed  of  sandstone,  shale  and  siltstone  from  which 
carbonaceous  material,  cross-stratif ication  and  ripple  marks 
are  conspicuously  absent.  The  Jurassic-Cretaceous  boundary 
coincides  approximately  with  the  contact  between  the  marine 
and  non-marine  sequences.  Fossils  were  collected  at  one 
locality:  a  Jurassic  pentacrinus  on  Packrat  Creek,  183  m 
below  the  Caaomin  Formation.  In  addition,  near  the  top  of 
the  Nikanassin  Formation  abundant  but  unidentifiable  plant 


remains  were  observed. 
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The  Nikanassin  is  unconformably  overlain  by  the  Cadomin 
Formation  (Irish,  1965,  p.  54) .  Using  a  method  similar  to 
that  described  by  Cruden  and  Charlesworth  (1966),  bedding 
orientations  from  both  sides  of  the  unconformity  were 
collected  at  five  localities  in  an  attempt  to  determine  its 
angularity  and  the  dip-direction  of  Nikanassin  strata  in 
Cadomin  time.  At  one  locality  there  was  no  significant 
angularity.  Apparent  angularities  of  28°,  22°,  9°  and  8° 
were  observed  at  the  other  localities  with  inferred  dip- 
directions  for  Nikanassin  strata  in  Cadomin  time  of  293°, 
321°,  72°  and  90°,  respectively.  The  absence  of  angularity 
at  one  locality,  the  large  apparent  angularities  at  two 
localities,  and  the  variation  in  inferred  dip-direction 
suggest  that  in  fact  the  sub-Cadomin  unconformity  is  not 
angular  and  that  the  apparent  angularities  observed  at  four 
localities  suggest  an  irregular  erosional  surface. 


Cadomin  Formation 


The  Lower  Cretaceous  Cadomin  Formation,  a  prominent 
ridge- forming  unit  whose  thickness  varies  from  22  to  38  m, 
is  made  up  of  a  lower  conglomerate  (8  to  15  m) ,  a  middle 
sandstone  (3  to  6  m)  and  an  upper  conglomerate  (9  to  18  m)  . 
The  conglomerates  consist  of  well-rounded  chert  and 
quartzite  pebbles  and  cobbles  in  a  medium-  to  coarse-grained 
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lithic  sandstone  matrix.  White,  grey,  black,  green,  banded 
and  occaisionaliy  red  chert  makes  up  approximately  60%  of 
the  rock.  The  quartzite  fragments  are  generally  white,  buff 
and  pinkish  in  colour.  The  lower  conglomerate  tends  to  be 
coarser  grained  and  to  have  more  green  chert  particles.  The 
conglomerates  are  silicified  and  extremely  well  indurated: 
fractures  generally  cut  across  fragments  rather  than  around 
them.  The  medium-grained  lithic  sandstone  contains  some 
fine-grained  silty  and  coaly  zones.  Poorly  preserved  plant 
remains  were  observed  at  the  base  of  the  formation  and  in 
the  sandstone  member.  The  upper  conglomerate  grades  into  the 
overlying  Luscar  Formation  within  a  zone  less  than  1  m 
thick. 


Luscar  Formation 


At  550  to  600  m  thick,  the  Luscar  Formation  is  the 
thickest  unit  present  and  underlies  more  than  half  the  study 
area.  Its  lithologies  are  many  and  varied  and  include 
sandstone,  shale,  siltstone,  conglomerate  and  coal.  Grey  and 
brown,  buff  weathering,  lithic  and  argillaceous  sandstone 
make  up  half  the  formation.  The  sandstones  in  the  upper  100 
m  are  grey  in  colour  and  are  thicker  bedded  than  those  lower 
in  the  formation.  Grey  tc  dark  grey  shale,  mudstone  and 
carbonaceous  shale  comprise  one  quarter  of  the  formation. 
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whereas  dark  grey  siltstone  accounts  for  nearly  20^. 
Approximately  2.5?c  of  the  formation  consists  of  coal  seams 
ranging  in  thickness  from  a  few  centimeters  to  4.5  m. 
Conglomerates  account  for  less  than  2%  of  the  formation. 
Concretionary  ironstone  hands  and  nodules  associated  with 
shale  and  siltstone  occur  in  minor  amounts.  Practically  all 
the  coal  occurs  within  the  lower  310  m  and  most  of  the 
conglomerates  within  the  lower  95  m  of  the  formation. 
Although  no  complete  section  through  the  formation  is 
exposed,  a  composite  section  based  primarily  on  exposures 
along  Packrat  Creek  and  Hells  Creek  has  been  established  and 
is  included  as  Appendix  8. 

four  map  units,  established  by  the  author,  are  present 
within  the  formation  (fig.  9)  and  are  described  below. 
Although  the  thicknesses  and  facies  of  the  horizons  within 
these  units  vary  somewhat,  in  general  they  are  continuous 
over  most  of  the  area.  Measurements  of  the  thicknesses  of 
the  units  at  different  localities,  although  too  few  to 
obtain  a  definite  rate  and  direction  of  thinning,  are 
consistent  with  a  northeasterly  direction  of  thinning.  The 
numbering  of  the  coal  seams  given  below  is  that  used  by 
McIntyre  Mines. 


Unit  A 

Coal  seam  #1,  1.2  to  1.8  m  thick,  occurs  20  m  above  the 
base  of  the  formation  and  is  directly  overlain  by  a 
resistant  4.6  m  thick  sandstone  and  co rglomerate  bed.  Coal 
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Stratigraphic  column  of  the  Luscar  Formation 


FIGURE  9 


4 


41 


W ,  ap  to  0.8  m  thick,  may  be  present  about  3  m  below  the  #1 
seam.  Coal  X/  also  less  than  1  m  thick,  may  be  present  about 
45  m  above  the  formation  base  and  is  overlain  by  a  7  m  thick 
sandstone  and  conglomerate  bed.  A  7  m  thick,  massive,  fine¬ 
grained  sandstone  is  located  60  m  above  the  base  of  the 
formation  and  3  m  above  this  the  0.5  m  thick  #1.5  coal  may 
be  present.  Located  approximately  85  m  above  the  base  is  the 
1  m  thick  #2  seam.  The  thicknesses  and  stratigraphic 
positions  of  all  except  the  #1  seam  were  obtained  from 
borehole  data. 


Unit  B 

Located  SI  m  above  the  base  is  a  31  m  thick  section  of 
marine  mudstone  and  shale  with  thin,  cross-bedded,  fine¬ 
grained  sandstones  and  rare  hard,  nodular,  silty,  orange- 
brown  weathering  mudstones.  This  unit  is  of  uniform 
thickness  and  is  persistent  over  the  entire  study  area. 


Unit  C 

The  most  distinct  marker  horizon  in  the  formation  is  a 
ridge- forming,  15  to  20  m  thick  sandstone  which  occurs  about 
150  m  above  the  base.  The  sandstone,  referred  to  locally  as 
the  *  Salt  and  Pepper  Sandstone1,  is  medium-grained,  lithic 
and  weathers  to  a  light  grey  with  some  reddish-brown 
limonite  staining.  Coal  seam  #3,  1  to  1.8  m  thick,  directly 


. 


. 


42 


overlies  the  sandstone  and  a  1.5  to  7.5  m  thick  pelecypod 
coquina,  referred  to  locally  as  the  'Clam  Zone',  begins  1.5 
m  above  the  coal.  Coal  seam  #4  is  located  nearly  200  m  above 
the  base.  It  ranges  in  thickness  from  about  4  m  on  the 
northwest  side  of  the  Smoky  River  to  between  3  and  3.3  m  on 
the  southeast  side,  and  reaches  a  maximum  thickness  of  over 
6  m  just  north  of  the  study  area.  A  massive,  fine-grained 
and  silty  sandstone  referred  to  locally  as  the  'Super  4 
Sandstone'  occurs  3  m  above  #4  seam.  Along  Two  Camp 
Anticline  this  sandstone  reaches  a  maximum  thickness  of  24 
m.  It  is  characterized  by  rapid  facies  changes,  and  at  some 
localities  it  is  absent  and  replaced  by  a  silty  and  shaly 
sequence.  Located  about  230  m  above  the  base  is  the 
discontinuous  0.2  m  thick  #5  seam  which  is  of  no  value  as  a 
marker  horizon  and  was  identified  in  only  a  few  boreholes. 
Coal  seams  #6  and  #7,  also  discontinuous  with  thicknesses  up 
to  1.2  m  and  separated  by  0  to  1.2  m  of  mudstone,  are 
located  about  240  m  above  the  base.  Seam  #8,  from  3  to  10.  7 
m  above  coal  #7,  has  a  thickness  of  1  to  1.4  m.  The 
stratigraphic  interval  between  #4  and  #8  seams  thins  from  50 
m  on  Mount  Hamell  to  35  m  on  Grande  Mountain.  Located  about 
295  m  above  the  base  is  seam  #10  (previously  named  #9  by 
Landes,  1962)  which  ranges  in  thickness  from  1  to  3  m. 

At  300  m  and  360  m  above  the  base  of  the  formation  are 
two  medium-  to  fine-grained  lithic  sandstones,  ranging  in 
thickness  from  15  to  24  m  and  12  to  21  m,  respectively.  The 
shaly  strata  separating  these  sandstones  are  by  no  means 
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uniform  in  thickness,  as  evidenced  by  borehole  6118  (26596E 
28105N),  located  on  Campbell  Flats,  where  the  shales  are 
absent  and  the  combined  thickness  of  the  sandstones  is  36  m. 
The  sandstones  weather  light  grey  with  large  oval,  reddish- 
brown,  limonite  stains.  Despite  being  only  poorly  indurated, 
the  sandstones  generally  form  prominent  ridges. 

Unit  D 

(Jnit  D  is  represented  by  upper  160  to  180  m  of  the 
formation.  A  granular  and  pebbly  sandstone  with  a  thickness 
of  18  to  25  m  is  located  about  400  m  above  the  base  and  160 
m  below  the  top  of  the  formation.  It  consists  of  a  lower 
granular  section  (7.6  m) ,  an  intervening  medium-grained 
sandstone  (5  to  10  m)  and  an  upper  granular  and  pebbly 
sandstone  (6.1  m) . 


Fossils 

With  the  exception  of  the  upper  120  m,  plant  fossils 
and  remains  are  abundant  throughout  the  formation.  In 
addition  to  plant  fossils,  the  following  important 
collections  of  pelecypods  and  gastropods  were  made. 

1.  Pelecypods  were  collected  from  the  coquina  horizon 
locally  referred  to  as  the  ■Clam  Zone'.  This  interval  is 
widely  distributed  and  can  be  found  overlying  the  #3 
coal  throughout  the  entire  study  area. 

2.  On  Grande  Mountain  (map  reference  37941E  1 446  IN), 


- 

. 


44 


very  veil  preserved  specimens  of  the  pelecypod  Unio  and 
some  gastropods  were  collected  along  a  road  cut  about  6 
m  above  the  #4  seam,  200  m  above  the  formation  base, 
from  the  silty  facies  equivalent  to  the  ‘Super  4 
Sandstone 1 . 

3.  On  Grande  Mountain  {map  reference  37826E  11037N), 
approximately  470  m  above  the  base,  some  gastropods  and 
peiecypods  were  collected  from  outcrops  along  a  cliff. 
This  fossil  zone  is  a  reliable  marker  horizon  within 
Unit  D  and  can  be  identified  in  borehole  samples  and  at 
other  outcrops. 


Regional  Stratigraphic  Correlation 
The  marine  mudstones  and  shales  located  about  90  m 
above  the  base  of  the  Luscar  Formation  probably  correlate 
with  the  Moosebar  Formation  of  northeastern  British 
Columbia.  There,  as  at  Grande  Cache,  the  beds  (belonging  to 
the  Gething  Formation)  that  underlie  the  marine  shales 
differ  from  the  overlying  strata  (belonging  to  the  Gates 
Formation)  which  have  thicker  and  more  abundant  coal  seams. 
At  Mount  Torrens,  65  km  to  the  northwest  of  Grande  Cache, 
Stott  (1968)  indicated  that  the  Moosebar  Formation  is  42  m 
thick  and  is  located  91  m  above  the  Cadomin  Formation.  At 
Grande  Cache  this  shale  wedge  is  significantly  thinner  than 
farther  north  and  continues  to  thin  southward  (Kilby,  1978). 
In  the  central  and  northern  foothills  of  Alberta,  thick 
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seams  observed  at  several  localities  occupy  similar 
stratigraphic  positions  within  the  Luscar  Formation.  Kilby 
(1978)  has  suggested  that  these  seams  may  be  correlated. 
Among  these  seams  are  the  #4  coal  seam  at  Grande  Cache,  the 
Jewel  Seam  at  Cadomin  and  Luscar,  the  Kennedy  Seam  at 
Mountain  Park  and  the  thick  seam  near  the  base  of  the  Gates 
Formation  at  Mount  Torrens  and  Mount  Belcourt,  100  km  to  the 
northwest  of  Grande  Cache  (Fig.  10). 

The  Mountain  Park  Formation  as  described  by  MacKay 
(1930)  is  not  distinguishable  from  the  Luscar  Formation  in 
this  area.  However,  the  uppermost  100  m  of  the  Luscar 
Formation,  composed  of  non-marine,  barren  sandstones  and 
siltstones,  is  possibly  equivalent  to  the  Mountain  Park 
Formation. 


Shaftesbury  Formation 


This  formation  has  limited  exposure  in  the  study  area 
and  is  mainly  restricted  to  the  Syncline  Hills  syncline.  It 
consists  of  140  m  of  dark  grey,  fissile,  marine  shale 
commonly  with  concretionary  ironstone  bands  and  some 
siltstone  and  sandstone  layers.  The  boundary  between  the 
Upper  and  Lower  Cretaceous  lies  near  the  top  of  the 
formation.  The  contacts  with  the  underlying  Luscar  and 
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FIGURE  10 


Comparison  of  sections  from  four  areas 
in  the  Rocky  Mountain  Foothills. 


overlying  Dunvegan  formations  are  sharp  and  gradational, 
respectively. 
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Dunvegan  Formation 


The  resistant,  ridge-forming  Dunvegan  Formation  is 
exposed  in  the  Syncline  Hills  syncline  and  consists  of  30  to 
45  m  of  interbedaed  sandstone,  shale  and  siltstone.  The 
basal  section  of  the  formation  contains  a  7.5  to  9  m  thick 
well  indurated,  massive  and  blocky  guartzose  sandstone  above 
which  lies  the  more  variable  interbedded  section.  The 
contact  between  this  formation  and  the  overlying  Kaskapau 
Formation  is  gradational. 


Kaskapau  Formation 


Exposures  of  the  Upper  Cretaceous  Kaskapau  Formation 
are  restricted  to  the  Syncline  Hills  syncline  and  to 
directly  below  the  Syncline  Hills  thrust.  Its  dark  grey, 
thin-bedded  and  fissile  marine  shales  commonly  are  silty  and 
contain  clay  ironstone  bands.  These  shales  are  very 
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incompetent  and  structures  involving  them  are  generally 
complicated  by  minor  faulting  and  folding,  making  thickness 
determinations  difficult.  However,  this  formation  is 
believed  to  attain  a  maximum  thickness  of  460  to  500  m 
(Irish,  1965) ,  of  which  only  the  lowest  100  m  are  present  in 
the  study  area. 


STRUCTURAL  GEOLOGY 


Strata  in  the  Rocky  Mountain  Foothills  of  the  Smoky 
River  region  are  folded  and  cut  by  numerous  southwesterly 
dipping  thrust  faults  (Fig.  1) .  Most  of  the  thesis  area  is 
underlain  by  Nikanassin  to  Kaskapau  strata  belonging  to  the 
Syncline  Hills  thrust- sheet.  To  the  southwest  these  strata 
are  bounded  by  Nikanassin  and  younger  strata  of  the  Cowlick 
thrust-sheet  and  to  the  northeast  by  Kaskapau  and  older 
strata  of  the  Mason  thrust-sheet.  The  area  has  been  divided 
into  22  domains  within  each  of  which  folding  can  be 
considered  cylindrical  (Fig.  11). 


Cowlick  Thrust-sheet 


Traceanle  for  about  100  km  along  strike,  the  Cowlick 
fault  is  one  of  the  larger  thrusts  in  the  Alberta  Foothills 
north  of  the  Athabasca  River  (Irish,  1S65).  Within  the 
Cowlick  thrust-sheet  only  the  Nikanassin  and  Cadomin 
Formations  between  Mount  Hamell  and  Grande  Mountain  were 
mapped.  The  fault  has  a  dip  of  approximately  60°  SW  beneath 
Mount  Hamell,  steepening  to  about  70°  SW  on  Grande  Mountain, 
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and  cuts  up  section  to  the  northeast  at  about  20°.  It 
superimposes  Nikanassin  and  Cadomin  on  Luscar  strata  and  has 
an  apparent  displacement  of  about  2.5  km.  Small  sub-parallel 
thrust  faults  with  displacements  less  than  10  m  are 
particularly  abundant  in  the  basal  50  m  of  the  thrust-sheet. 
Two  macroscopic  and  several  mesoscopic  chevron  folds  with 
axial  planes  dipping  steeply  northeast  and  hinges  plunging 
gently  northwest  were  observed  on  Mount  Hamell. 

On  Mount  Hamell  the  plunge  of  folds  in  the  thrust-sheet 
is  to  the  northwest.  In  the  floor  of  the  Smoky  River  the 
plunge  becomes  southeastwards  and  steepens  to  10° ,  bringing 
the  Cadomin  Formation  in  contact  with  Luscar  strata  across 
the  Cowlick  fault.  On  Grande  Mountain  the  plunge  is  again  to 
the  northwest. 


Syncline  Hills  Thrust-sheet 


The  northeastern  limit  of  the  study  area  is  marked  by 
the  Syncline  Hills  Thrust,  dipping  to  the  southwest  at 
approximately  50°.  The  fault,  which  brings  Luscar  strata 
onto  the  Kaskapau  Formation  in  the  footwall,  is  not  a  major 
thrust  and  cannot  be  traced  along  strike  for  any  great 
distance . 

Most  of  the  Synciine  Hills  thrust-sheet  is  underlain  by 
the  Campoell  Flats  anticlinorium  (Plate  1).  The 
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PLATE  1.  Southeastern  slope  of  Mount  Hamell  viewed  from  Grande  Mountain;  N-Ni kanassin ; 

Ca-Cadomin;  A,  B,  C3  D-Luscar  units  A,  B,  C,  D  respectively;  F-Cowlick  Thrust 
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anticlinorium  exposes  Nikanassin  and  Cadomin  strata  along 
the  Smoky  River,  with  Luscar  strata  forming  the  higher 
ground  on  Mount  Hamell  and  Grande  Mountain.  The  Shaftesbury, 
Dunvegan  and  Kaskapau  Formations  are  confined  to  the  higher 
ground  on  both  sides  of  the  river,  and  to  the  core  of  the 
Syncline  Hills  syncline  which  borders  the  Campbell  Flats 
anticlinorium  to  the  northeast  (Fig-  12). 

At  the  level  of  the  Nikanassin  and  Cadomin  Formations 
the  Campbell  Flats  anticlinorium  resembles  a  box-fold  with  a 
horizontal  central  limb.  The  Cadomin  Formation  in  the  fold 
which  borders  the  central  limb  to  the  southwest,  known  as 
the  Sterne  Creek  anticline,  actually  projects  about  250  m 
above  the  level  of  the  Cadomin  in  the  central  limb,  forming 
an  ear  to  the  box-fold-  Its  northwesterly  trending  fold-axis 
is  horizontal  in  domains  3  and  4  and  gently  plunging  to  the 
northwest  in  domain  5-  There  is  a  corresponding  change  in 
the  orientation  of  the  axial  plane  from  37°  54°  to  46°  32° 
and  in  the  value  of  the  apical  angle  from  95°  to  60°  (axial 
plane  orientations  given  as  dip  direction  and  dip,  fold-axis 
orientations  given  as  trend  and  plunge) . 

Two  Camp  anticline,  the  fold  bordering  the  central  limb 
of  the  anticlinorium  to  the  northeast,  is  a  simple  chevron 
fold  at  the  level  of  the  Cadomin  Formation,  in  domains  9  and 
10-  Towards  the  northwest  in  domain  8,  within  the  Luscar 
Formation,  the  fold  becomes  more  complex,  forming  a 
northeastern  ear  to  the  Campbell  Flats  anticlinorium-  The 
southeasterly  trending  fold-axis  is  essentially  horizontal 
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in  domains  9  and  10  nut  steepens  to  15°  SE  in  domain  8, 
whereas  the  dip  of  the  northwesterly  striking  axial  surface 
changes  from  60°  SW  in  domains  9  and  10  to  75°  SW  in  domain 
8.  There  is  a  corresponding  decrease  in  the  apical  angle, 
from  120°  to  80°. 

Fox  Creek  syncline  and  a  small  unnamed  anticline  lie  to 
the  northeast  of  Two  Camp  anticline.  The  fold-axes  of  these 
rounded  chevron  folds  plunge  to  the  southeast  at  5°  and  the 
axial  planes  dip  to  the  southwest  at  75°.  Along  strike  to 
the  northwest  the  two  folds  are  replaced  by  a  steep, 
northeasterly  dipping  reverse  fault  (Fig.  13) .  In  the 
hanging  wall  of  the  fault,  the  tight,  chevron.  Syncline 
Hills  syncline,  with  an  apical  angle  of  50°,  plunges  gently 
to  the  west-northwest  with  an  almost  vertical  axial  plane. 
Directly  below  the  fold  to  the  northeast  is  the  Syncline 
Hills  thrust. 

The  major  box-fold  so  conspicuous  on  Mount  Hamell, 
cannot  be  identified  southeast  of  the  Smoky  River  on  Grande 
Mountain.  However,  Sterne  Creek  anticline  and  Syncline  Hills 
syncline  are  easily  located. 

Sterne  Creek  anticline  with  a  horizontal,  southeasterly 
trending  fold-axis  has  an  axial  plane  orientation  of  50°  78° 
and  an  apical  angle  of  about  100°.  The  fold  is  no  longer  a 
simple  structure,  however,  and  in  its  southwestern  limb 
contains  both  southwesterly  and  shallow  northeasterly 
dipping  thrusts  which  cut  up  section  at  low  angles  to  the 
northeast  and  southwest,  respectively.  Southeast  of  the 
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crest  of  Grande  Mountain,  Sterne  Creek  anticline  appears  as 
an  asymmetric  chevron  fold,  disrupted  on  the  southwest  limb 
by  a  small  thrust  with  a  displacement  of  about  1  km  that 
brings  Cadomin  and  Nikanassin  beds  onto  Luscar  strata. 

Between  Sterne  Creek  anticline  and  the  continuation  of 
Fox  Creek  syncline,  four  small  folds,  which  die  out  on  the 
lower  slopes  of  Grande  Mountain,  were  observed.  Two  Camp 
anticline  cannot  be  correlated  across  the  river.  Fox  Creek 
synciine  and  the  adjacent  anticline  also  die  out  below  the 
crest  of  Grande  Mountain.  All  these  folds  have  nearly 
horizontal  southeasterly  trending  fold-axes  and  axial 
surfaces  that  dip  southwest  at  about  70°.  Syncline  Hills 
syncline  continues  across  the  northern  slopes  of  Grande 
Mountain  as  a  more  open  fold  with  an  apical  angle  of 
approximately  100  .  The  geographic  extent  and  geometry  of 
the  folds  on  individual  surfaces  are  displayed  on  structure 
contour  and  overburden  maps  produced  for  coal  seams  3,  4  and 
10  (Fig.  14,15,16,17,18,19). 

The  structure  of  the  area  southeast  of  Grande  Mountain 
remains  obscure  because  of  poor  exposure  and  the  virtual 
absence  of  drill  hole  data.  A  quarry,  beside  highway  40,  in 
the  faulted  and  folded  strata  of  the  Cadomin  Formation 
provided  good  artificial  outcrops.  The  geometry  of  the 
folding  is  that  of  a  box-fold  with  one  ear  on  the  northeast 
limb  which  is  faulted  by  two  northeasterly  dipping  thrusts 
with  a  combined  displacement  of  about  250  m.  The 
northwestern  limb  of  the  box-fold  is  a  continuation  of 
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Sterne  Creek  anticline,  with  the  ear  structure  not  evident. 


Mesoscop>ic  Structural  Geology 


Minor  bedding-plane  faults  with  almost  no  apparent 
displacement  are  common  in  the  interbedded  strata  of  the 
Luscar  and  Nikanassin  Formations.  The  faults  follow  bedding 
for  a  distance,  then  cut  up  section  for  generally  a  meter  or 
less  and  then  resume  a  position  parallel  to  another 
stratigraphic  level.  Larger  faults  with  displacements  of  1 
to  10  i  are  common  (dipping  both  to  the  southwest  and 
northeast)  and  normally  make  an  angle  with  bedding  of  30  or 
less. 

Most  of  the  mesoscopic  folding  is  congruent  with  the 
macroscopic  folds  although  incongruent  folds  do  occur 
adjacent  to  faults  (Fig.  6) .  Tight  chevron  folds  and 
numerous  kink  bands  varying  in  width  from  0.5  to  5  m  are 
common.  Insufficient  data  were  collected  on  jointing  and 
coal  cleat  to  warrant  comment. 


Discussion 


Abundant  bedding  plane  slickensides  perpendicular  to 
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the  fold-axes  throughout  much  of  the  area  indicate  that  the 
flexural  slip  mechanism  of  deformation  was  active.  This 
mechanism,  along  with  the  narrow  hinge  zones,  planar  limbs 
and  the  discontinuity  of  the  folds  both  along  strike  and 
vertically,  suggest  that  kink  folding  on  a  macroscopic  scale 
was  the  primary  means  of  deformation  (Faill,  1969,  1973). 

The  sudden  changes  in  orientation  of  many  fold-axes  and  the 
limited  geographic  extent  of  the  cylindrical  domains  (e.g. 
Two  Camp  anticline)  are  characteristic  of  kink  folds. 
Numerous  mesoscopic  kink  bands  were  also  observed.  The 
geometry  of  most  of  the  folds  in  the  area  can  be  described 
in  terms  of  large  scale  kink  bands  of  various  sizes  and 
orientations  (Fig.  20) . 

The  steep  Cowlick  and  Syncline  Hills  faults  are 
probably  rotated  imbricate  thrusts  in  the  hanging  wall  of  a 
major  thrust,  possibly  the  Mason  thrust,  which  is  the 
largest  fault  in  the  foothills  north  of  the  Athabasca  River. 
This  fault  could  also  be  acting  as  a  zone  of  decoilement 
between  the  Jurassic  and  the  pre-Jurassic  strata. 

Three  levels  cf  disharmonic  mesoscopic  deformation  are 
observable  in  the  stratigraphic  succession.  Much  of  the 
deformation  present  in  the  Nikanassin  Formation  terminates 
near  the  base  of  the  Cadomin  Formation.  The  second  level  is 
bounded  above  by  the  thick  shale  horizon  located  90  m  above 
the  base  of  the  Luscar  Formation.  The  third  level,  within 
which  minor  faulting  predominates  over  folding,  lies  in  the 
upper  Luscar  strata  above  the  shale  horizon.  From  the 
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FIGURE  20 
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geological  map  (Fig.  4) ,  many  folds  and  faults  can  be  seen 
to  terminate  or  alter  in  configuration  near  these 
boundaries. 


SUMMARY 


This  study  has  attempted  to  contribute  to  the 
understanding  of  the  geology  of  the  foothills  of  west- 
central  Alberta  aDd  to  provide  some  insight  into  the 
applications  of  some  computer-based  techniques  that  may  be 
used  during  exploration  for  coal  seams  in  structurally 
complex  areas.  Specifically,  this  thesis  accomplishes  the 
aims  outlined  in  Chapter  1,  namely:  (1)  to  determine  the 
stratigraphy  of  the  Luscar  Formation;  (2)  to  describe  the 
structure  of  the  Luscar  and  adjacent  formations  in  the 
vicinity  of  Mount  Hamell  and  Grande  Mountain;  (3)  to  develop 
a  computer-based  structure  contouring  procedure  that 
utilizes  orientation,  positional  and  stratigraphic  data;  and 
(4)  to  d  emonstrate  the  usefulness  of  some  computer-based 
techniques  in  determining  and  analysing  the  structural 
geology  of  coal  measures.  The  structure  contouring  procedure 
is  particularly  useful  during  preliminary  exploration  for 
coal  to  assess  potential  mining  areas  and  to  prepare 
structure  contour  and  overburden  maps  of  coal  seams  without 
the  need  for  extensive  drilling. 

Additional  detailed  studies  such  as  this,  through 
cooperation  between  the  university,  industry  and  government 
research  agencies,  are  required  to  gain  a  fuller 
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understanding  of  the  geology  of  the  Foothills  of  the  Western 
Canadian  Rocky  Mountains  and  to  develop  new  and  refine 
existing  computer-based  procedures  helpful  in  determining 
and  analysing  their  structure. 
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APPENDIX  1. 

The  format  and  contents  of  the  file  DATA  are  as 
follows: 

LINE  1 

Column  Description 


1  Numeric  code  specifying  form  of  orientation  data: 

1=  trend  and  plunge 
2=  dip-direction  and  dip 
3=  strike  (right-handed)  and  dip 


2-41 

Title . 

SECOND 

TO  SECOND  LAST  LINES 

Column 

Description 

1-4 

Identif ication 

number. 

9-10 

Binary  code  of 

which  the  first  digit  identifies 

formation  and  the  second  stratigraphic  position 
within  the  formation: 

90=  Kaskapau  Formation 

80=  Dunvegan  Formation 

70=  Shaftesbury  Formation 

40=  Luscar  Formation  (undivided) 

57=  upper  100-125  m  of  the  formation 
56=  25  m  interhedded  medium  grained  and  pebbly 

sandstones 

55=  35  m  of  sandy  and  shaly  strata 

54=  18  m  lithic  sandstone  (upper  sandstone) 

53=  40  m  of  interbedded  shale  and  sandstone 

52=  20  m  lithic  sandstone  (lower  sandstone) 

51=  Coal  seams  10,  8,  7,  6  and  related  strata. 
50=  15  m  of  shale 

49=  Coal  seam  4,  "Super  4"  sandstone  and 
related  strata 

48=  15  m  of  shale  and  sandstone 

47=  Coal  seam  3,  "Salt  and  Pepper"  sandstone 
and  related  strata 

46=  25  m  of  mudstone,  siltstone  and  sandstone. 

45=  31  m  of  mudstone 

44=  Coal  seams  2,  1.5  and  related  strata. 

43=  Coal  seam  X  and  related  strata. 

42=  Coal  seam  1  and  overlying  sandstone 
41=  lower  20  m  of  the  formation  with  coal  W 


68 


' 


-  . 


69 


30=  Cadomin  Formation  (undivided) 

33=  upper  conglomerate 

32=  middle  sandstone 

31=  lower  conglomerate 

20=  Nikanassin  Formation 

11-14 

Distance  in  feet  of  exposed  horizon  above  base  of 
formation  to  which  it  belongs. 

15-29 

X  (easterly),  Y  (northerly)  and  Z  (elevation) 
coordinates  in  feet. 

30 

Station  position  reliability  indicator  found  by 
dividing  estimated  possible  error  in  feet  by  10- 

31-38 

Presence  of  additional  information  recorded 
elsewhere  was  indicated  by  a  1=yes  or  0=no  entry 
in  these  columns. 

31=  bedding  plane  linears 

32=  mesoscopic  folds 

33=  faults 

34=  joints 

35=  lithology 

36=  fossils 

37=  general  observations  including  sketches 

38=  photos 

39-79 

Bedding  orientations  recorded  as  six  6-digit 
entries  of  dip-direction  and  dip  with  a  space 
between  each  entry-  Dip-direction,  the  azimuth  of 
the  direction  normal  to  the  strike  in  which  the 
strata  young,  can  range  from  0  to  400  grads-  Dip, 
the  angle  through  which  the  strata  have 
apparently  been  rotated,  can  range  from  0  to  200 
grads,  being  greater  than  100  grads  in  the  case 
of  overturned  beds. 

80 

Presence  of  additional  bedding  orientations 
recorded  elsewhere  was  indicated  by  a  1=yes  or 
Q=no  entry  in  this  column- 

LAST  LINE 


Column 

Description 

1  * 

/ 

. 


f*  -  ;  c  t  4xt 


■ 
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APPENDIX  2. 

The  format  and  contents  of  the  file  BOREHOLE  are  as 
follows: 

FIRST  TO  LAST  LINES 


Column 


Description 


1  Alphanumeric  code  designating: 

5=  topographic  surface 
1-4=  coal  seams  1-4  (tops) 

6-9=  coal  seams  6-8  and  10  (tops) 

D=  Cadomin  Formation  (top) 

Z=  fault 

2-16  Coordinates  of  intersection  of  above  surface  and 
borehole  as  in  DATA. 

71-74  Identification  of  which  digits  1  and  2  indicate 

year  of  drilling  and  3  and  4  identify  borehole  in 
McIntyre  Mines*  records. 


80 


1  (see  program  DRAH) . 


. 
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APPENDIX  3. 

Listings  of  computer  programs  written  for  data 
processing  with  appropriate  inputs  and  outputs. 

DELET 


input  5=DAf A 

output  6=-DAI A 1 


1 

C 

— 

TO  COM  VI.1.T  GliAOS  T1 

2 

c 

— 

‘JO  STRATIGRAPHIC  I'.i 

3 

L>  IM  i'W  TIGS  X  (50) 

4 

INTKGiif  m,Db,X 

5 

2  Si  T hC£r,  CO  J  L  ( 26 )  /  '  A 

6 

* • L ' , »  tt ' , ' V  * , • 0 * ,  *  P  * 

.  7 

K  =  0 

fi 

c 

— 

HEAD  1U  Si  AW  DATA 

n 

1 

F.  E  Ali  ( 5 , 1 0  ,  hH  L>=3 )  1*0, 

10 

*  (X  (J)  ,2  =  1 ,  12)  ,HO 

n 

c 

— 

CliAWGL  UliklJS  To  uF' 

12 

UJ  2  J=1 , 12 

13 

r.  (j)  <j)  *o  .9*o  .5 

14 

COFTI.N'JL 

15 

c 

— 

ASSIG.V  TiiL  LET  TEH  < 

16 

If  [it*  .EC.  0)  Go  TO 

17 

IF  (fa  .40.  K)  Gu  TO 

18 

i 

19 

.*-fa 

20 

c 

- - 

»  FITE  THE  Y.L4  I/XT  A 

21 

C 

\'R  IT  E  (n  ,  11)  COi/i  (5 )  , 

22 

(X  (J)  ,J*1,  12)  U 

23 

GG  TO  7 

24 

* 

o 

n  J*  IT T  (6 , 1  /. )  :I o ,  F t  ,  L)ii 

25 

-  (A  (J)  ,3  =  1  ,  12)  ,ai, 

2b 

7 

CGXTI  ;«U  E 

27 

GO  TO  1 

20 

1C 

f Ortfl/T <14 ,41,12,14, 

29 

1 1 

f  0  HR  AT  (  A  1  ,14, 3X  ,12, 

30 

12 

FOHEAT ( IX ,14,32,12, 

31 

3 

STO?- 

32 

ESI/ 

SD  OF 

FILE 

l  c;  i 


»  »  «  n «  i  t'i  i  t*'  '  vi  »  y  •  171/ 

*»*•<•»  t  i  )  it  / 


. 


SORT 


input 

5=-DAI A  1 

output 

6=-DATA2 

1  C 

— 

TO  SORT  DATA  kS  IN PDT  FOR  EEAK 

2 

DIMENSION  X  (50) 

3 

INTEGER  A,B,C,X 

4 

WHITE  (6,30) 

5 

DO  1  1=1,800 

6 

R2AD  (5, 10  ,END=3)  K,  A,  5,0,  (X  (J) 

7 

WRITE  (6,20)  K,A,3,C,  (X(J)  ,J=1, 

S 

1 

CONTINUE 

Q 

> 

10 

FORMAT  ( A  1 , 131,315, 9X, 213,5 (IX, 

10 

20 

r‘0 Ril  A T  (Al  ,3 15 , 6  ( 1 1 , 2X3 )  ) 

11 

30 

POSCAT( *2* ,3X, «GRANDE  CACHE1) 

12 

40 

PORK  A  T  (»/’) 

13 

3 

CONTINUE 

14 

WRITE  (6,40) 

15 

STOP 

16 

END 

2ND  OP  PILE 

■ 
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HSS08T 


inputs  4=-DAT A  1 

5=-*MEAN  (from  MEANS) 
7=-KV  (from  MEANS) 

output  6-MEANN 


1  C - TO  ARRANGE  DATA  AS  INPUT  INTO  HUflERICAL  PROGRAMS 


2 

INTEGER 

,U,F,S, 

? 

3 

'A  RITE 

« 

4 

DO  1  I»1,iSuO 

5 

■  F.LAD  (5, 20  ,  ESD* 

■-2)  K,X 

,i, 

,2, 

L>,  P 

b 

HE. 10(4, 30)  K  ,  i 

p  1  ,*♦  \  y 

1  ,f 

12, 

tt  3 , 1  4  ,  (-  5  ,  fl  6  ,  f.  7  ,  F.  H  ,  .1  '•< ,  **  1 0 

1 

HEAD  (7,60)  r.s , 

rr  1  1 

r  W  V-  t  i*  1  / 

L2 

XKITE  (t ,40)  K, 

t  ^  9  Y  9  i*  t 

0 ,  V ,  f» 

?!  ,C.:,L  1,L2, 110,  r«1,  .12, 13, 

9 

*  %  4  ,  fib  ,  P  6  ,  Y.  7  ,E3  ,  E  9  ,  If 

K  ,  1 

•  9  * 

10 

1  COSTIRD  h 

n 

10 

r’OEflAT('2'  ,3X,  •GRAND;';  CACUD') 

12 

20 

fOSiiAT  (A1 ,315,  IX  ,  21  3) 

13 

20 

?C  Ril  A  T  (14 , 4  x ,  1 2 , 1 1 ,  1 5  X  ,  9  11  ,4  11,11) 

14 

40 

FOR.iAT  (A  1 ,  Jl 5, 12 , 21  3, 2X,  11,  2X,  5 a  .2, 21, 213, ,  i  1 , 

15  • 

*2:c ,  1 1 , 2x , j  j;l  ,2::  ,i 4 , 4 x  ,14 , 2x  ,12) 

16 

60 

FORMAT (26X,I1 ,2X ,fb.2, IX ,13, U,x3) 

17 

50 

KOK;1AT  { V*  / 

18 

2 

CONTINUE 

19 

7SITF (6 ,50) 

20 

STOP 

21 

END  0 ?  FILS 

END 
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SELECT 


inputs  4=HEANN 

5=0UTNUM  (outcrop  numbers  of  data  stations  to  be 
retrieved) 

output  6=D0MX  (domain  file) 


1 

2 
3 

A 

5 

£ 

7 

9 

o 

10 
1  1 
1  2 

13 

1  A 

15 
1  6 
17 

19 

lo 

70 
21 
22 
23 
?  a. 

25 

26 
2^ 
29 

20 

gN  P  HP 


TO  SELECT  OATA  STATIONS  Fq0„  th£  fi,_e  MEAnh  . 

ANO  TO  PLACE  them  in  a  rlL5  D0MX  F09  OOMMN  . 

ANALYSIS  .  THE  STATION  NUMOEPS  APS  INPUT  INTn 
SELECT  IN  A  FILE  -OQMX. 

OIMFNSI3N  N(1 SO ) 


Z  t.oo  j  oj_iU  •  l  s  ,  rj , 


pFAL  C° 

OF  AO ( 5 •  1  3 )  9,C 
WPIT£{6. 131  ra,c 
00  1  J=2,150 


xPJJL 


gu.F 


S!FA0(5,10.EN0-=2)N(J) 


-$—■  ZA  l-»  ENa=?.JrLf  XjY_»  Z  »0U  ,  p  T  ,  nn  1  cc  .  L  .  L  1  ,  '»  1  .  M  ?  ,  »*3  .  MA  .  Mg.  . 


*M6  ,  M7  ,MJJ  ,  MP  ,M10  ,0,C  ,00 
IF (  N(  J)-C> A  ,6. A 

6  WPTT£(6,11  )<,X,Y.  !,OOi'l[,f|NtCc,i.,i.j,Mi,M?,M3,<H,w5, 
*M ,  M7  ,  M9  ,  MO  ,N1  0  t  9  tC  ,0>J 
GO  TO  1 

A  CON T I NUF _ 


GO  TO  5 

1  CONTINUE 

2  CONTINUE 
WPITEl 6. 15) 

10  FOPMAT(IA) 


11  PQPyATlAL,J3.LS.J.l  X  ,2M.  2/.  ,  £  j  f 


*2X  » 81 1 «2X* IAt AX i  I  a,  2X » f  7  j 
13  FOPMATI  I  1  .7H00MA IN,  13 ) 

15  FQQMAT ( •/• ) 


2X-»£a«2»2X  «?TJL . ?X ,  ? 1  ,  2X  .  I  1  . 


STOP 

ENO 


■ 


SLICKS 


input  5=PITCH  (file  containing  orientation  cf  planar 
surface  and  pitches  of  slickenside  striae) 

output  6=SLICKENS 


1 

C - 

SLICXS  IS  ul  SIU.JSD  T 

2 

C - 

LI HEAfvS  K&ASiJn£D  AS 

3 

INTKUSi?  DD,DI,PIT,MO, 

4 

DTK=1 .7543292-^ 

5 

3?{D-‘iO.O-»OTH 

6 

ftrUD  (5  ,  10)  X,  I  ,Z 

7 

2r:»0  (5,11)  ylj,  or 

6 

ifilTS  (6 ,12)  j;o,X,  It  ,L 

9 

.4  ; IC5  (o,13) 

10 

X'.IU =L>i)  ■*!)'?* 

11 

12 

no  i  i-i,5o 

13 

&SAu  (5,  i« ,  :::ib=2)  pit 

14 

X:»IT=PIT*0TL 

15 

xp=asi;;  (Cos  (r.  iu-xpn ) 

16 

?.C=il sD-ASI.M  (TA  J  (1£P)*T 

17 

It  (PIT-50)  4,5,5 

Ifi 

4 

:ct  -xou-xc 

19 

GO  10  4 

20 

5 

KT=XL*U-*-XC 

21 

6 

CuSTi..Mit: 

22 

IHii:*XY/l*Th*v  .3 

23 

IF  (Tfi  t'N  .LT.  <)) 

24 

I  F  (Th K;\’  .GI.  360)  TLEJi 

25 

PLDN-XI /DT  -^O  .5 

26 

:i h IT f  ( o  ,  1 1 )  Th  FA  ,  r-  LU :i 

27 

1 

CG!»TIillf2 

26 

4a 

CONTINUE 

29. 

9 al'I F.  (6,15 1 

30 

10 

Fvl  t  AT  («;  a,  13,  1a  ,  3  15) 

31 

1 1 

FO  i: «  A  T  (213) 

32 

12 

r  uhtik’l  CSTATIor:  •  ,  IX,  1 

33 

13 

r-Ot  >*AT  l  •  Til  5-1  • ) 

34 

14 

FOpOAC  (1  j) 

o5 

15 

FOLftAT { ‘LSD  •) 

3  6 

STOP 

37 

F‘;b 

PJ.'iC  TRi.NO  AMD  'J  ’ 

CLOCK  <*  IS E  PITCS  0;j  •*, 


.1 


3  oO 
3  6C 


2ND  OF  PILE 


. 


' 


. 


. 

APPENDIX  4. 

The  format  and  contents  of  the  file  MEANN  are  as 
follows: 

LINE  1 

Same  as  for  DATA. 

SECOND  TO  SECOND  LAST  LINES 


Column 

Description 

1 

Alphanumeric  code  designating  stratigraphic 
position,  coal  seam  number,  etc: 

Y=  Kaskapau  Formation 

X=  Dunvegan  Formation 

W=  Shaftesbury  Formation 

E=  Luscar  Formation,  interval  unknown 

F-V=  Luscar  Formation,  intervals  41-57 

1-4=  Luscar  Formation,  coal  seams  1-4 

6-9=  Luscar  Formation,  coal  seams  6-8  and  10 
B— D=  Cadomin  Formation,  intervals  31-33 

A=  Nikanassin  Formation 

Z=  Fault 

5=  surface  coordinates  of  a  borehole 

2-16 

Station  coordinates  as  in  DATA. 

18-23 

Mean  bedding  orientation. 

26 

Number  of  repetitive  bedding  orientations. 

29-36 

Concentration  parameter. 

39-44 

First  bedding  orientation  measured  at  outcrop. 

47 

Same  as  column  80  in  DATA. 

50 

Same  as  column  30  in  DATA. 

53-60 

Same  as  columns  31-38  in  DATA. 

63-66 

Same  as  columns  11-14  in  DATA. 

71-74 

Same  as  columns  1-4  in  DATA. 

77-78 

Same  as  columns  9-10  in  DATA. 

LAST  LINE 

Same  as  last  line  in  DATA. 


. 


:<.t. 


■ 


77 


APPENDIX  5. 

Listing  of  the  program  DRAW  with  the  appropriate  run 
commands- 

inputs  4=*S0URCE* 

xmin,  ymin  =  coordinates  of  lower  left-hand 
corner  of  map 

xrnax,  ymax  =  coordinates  of  upper  right-hand 
corner  of  map 

scale  =  scale  in  feet  per  inch  on  map 
ht  =  height  of  characters  in  tenths  of  an  inch 
dtic  =  distance  in  inches  between  distance 
markers  on  the  map  boundaries 

ocnum  =  0,  outcrop  numbers  not  plotted,  or  1, 
outcrop  numbers  plotted 
5=MAPPILE  (input  file  to  be  plotted) 

output  9=plotted  map 


-r-  ^ 


■ 


1 

2 

3 

4 

5 

6 

7 

a 

9 

70 

71 

1-t 

13 

14 

15 

1o 

17 

IB 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

3o 

37 

38 

35 

40 

41 

42 

4 3 

44 

45 

4b 

4  7 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

55 

oO 


78 


C***»****DliA«  IS  DESIGNED  TO  PLOT  DIP—  0  la  ECHO  K  AND  DIP 
C********S'll‘.iQlE  OK  A  tip  xLOHG  VIT3  BOUNDARIES. 

IKTEGL4  X  ,1  ,Xfl,  23  ,ID£KT 

REAL  ET,X3,2.»,Db,DI,N  ,BX  ,  bX  ,  OCHIJB 

CALL  PLOTS 

CALL  XLIEIT  (110  .0) 

32  AD  (4,  11)  XEIH,  2SIK  ,  XAAX  ,  XflA  X, SCALE  ,  2T,D  TIC „GCH0S 

XLLH  =  (XEAX—X2IN  )  /SCALE 

1  LLK=  { I  MX -XX IH) /SCALE 

IBOl»KD=XLB£*1  .0 

X5GUKD=T  LEK+G.5 

It  IK  1=X2Ia/1GQ 

yfiiNi=i  ciK/ioo 

SCALE!— SCALE/100 

C»*****«*  DA An  fcOOHDAXlES  ***** 

1000  CALL  AX2ES- (2.0, 1,-1) 

CALL  AXIS2  (1.0,  0.5,  *  *  1 ,  ILLS  ,  0 . 0  ,  I-INl  ,  SCALE!  ,DTIC) 

CALL  AXIS 2  (X30QKD ,0.5,'  •  1 ,2  LEN ,90.0 ,7  MM 1 rSCAL31 , DTXC) 

CALL  AIIS2 (1.0, 0.5, *  * , 1 , ILEA , 90 . 0 , in  la  1 , SCALE! , DTIC) 

CALL  AXIS2  (1 . 0  ,  IEOUXD  ,  1  *  ,  1  ,XLLN  ,0.0  ,1S.I  S 1  ,  SCALE!  ,DTI  C) 

C********  HEAD  OBIEfi'1  aTIOHS  AiiD  PLOT  SXEBOLS  ***** 

5  BEAD  (5, 10, END— 2)  X,T,bD,Dl,H ,IDSNT 
23=1 
13=1 

X3=  (X3-  (XKIN-SCALE)  )  /SCALE 
X  3=  (X  3-  (X  Eli* -SCALE  *0  .  5)  )  /SCALE 
1/ <X3  .LE.  1.)  GO  TO  1 
17  (X3  -GB.  X300ND)  GO  TO  1 
I?  (23  .LE.  .5)  GO  TO  1 
lr  (13  «G£.  XBGOhD)  GO  TO  1 

C*******«rt;2TEEElK£  I.P  STATION  TO  BE  PLOTTED  IS  AM  OGTCBOE, 
C*****»**j3oRESOLE, FOLD-AXIS  03  SLICXEK  SIDE 
17  (IDEKT-1)  300,400,500 
300  CALL  DBA* (13,23, X,ZT, DO, Dl,OCHDa) 

GO  TO  51 

400  CALL  S12rGL(X3,Y3,HT,010,u.G,-1) 

GO  TO  51 

500  CALL  CESO  (X3 , 73 , DD, DI,ET , IDE3T) 

51  CGLIIKUE 

1  COKTIKUE 
GO  1*0  5 

2  CALL  PLOT  (G.,0.,  999) 

10  70KAAT (IX ,15,15, oX, 73-0,73. 0 ,47X, 74. 0,5X, II) 

11  ?Ck£  AT  (87b  .  1) 

STOP 

EH  D 

£****•»* *********  SUBKOGTIjiL  DBAS  *■********»»*'»-*« 

.  SOfattOOl'lDE  DBAS  (X  ,  2 , 3  ,BI ,  DD,  DI  ,GC*  U2) 

SEAL  H 

BEAL  DD,DI,DDD,Z,7,» 

17  (DI  .EQ.  0.0)  GO  TO  50 
BEAL  I,I,S,ST,K,£It,P,Q,T 

C** ******  CHECK.  17  BEDS  OVrBTUSKED,  I?  SO  ADD  160  TO  DIP-DIEECTIOft 
IF  (DI  .LE.  90.)  GO  TO  15 
DDD=DDf130 . 

17  (ODD  .LE.  360.)  Go  TO  3 
DDD=DDb— 360 . 

V=A8S  (DCD-4  50.) 

GO  TO  6 

3  W=AES  (DDD-450.) 
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6  G=W+9G. 
ii=«f27U. 

K=2T/4. 

LD=LLD 

C»=******  COk'VSBT  DEGkLES  lu  RADIANS  **** 

15  Z=DD*1.74S329E-2 

U  =  (DD-t-lbu.)  *1  .745329E-2 
?=  (LD+90.)  *1.7453291-2 
KMDD-yO.)  *1.7453296-2 
I?  (DI  -LZ.  9u . )  GO  TO  20 

C********  COORDINATE  CALCS.  EOa  OVEHTuKfiZD  CIRCLE  SIE30L 
C********  A A’D  ITS  PLOTTING 
Xl=X+SIi.  (0)  *  (ST*  .33) 

1 1  — Y  +  COS  (D)  *(ST*.33) 

X2=X-*SIN  (V)  *  (ET*  .5)  • 

12=1  +  COS  (V)  *(ST*.5) 

X3=X2+S1N  (U)  *  (SI*  .33) 

I3=Y2+CQS  (0)  *  (HI*. 33) 

CALL  PLOT  ( X  ,  Y ,3) 

CALL  PLOT  (X  1,  Y  1 , 2) 

CALL  CIRCLE  (XI,  Y1  ,G,H,R,E,O.Q) 

CALL  PLOT  (X3,T3,3) 

CALL  PLOT  (J2,X2,2) 

DI=180.-D1 

C********  DETERSIVE  COORDINATES  EGG  S YE SOL,  LIP  AND 
C*****»»*  OUTChOP  NUSBER 
20  E=X*SIN  (Z)  *  (HT*.66) 

E=1+CQS  (2)  *  (3T*  .66) 

A=X+SIH  (Y)  *uT 
E  =  Y+COS  (7)  *HT 
C=X+SI5  (K)  *3T 
D=Y+COS  (K)  *ST 

IE  (2  .GT.  2 .35o 1942  .AND.  2  .LT.  3.9265903)  GO  TO  10 
IE  (2  .GE.  3.9269503  .AML.  2  .LT.  5.4977b6*«.)  GO  TO  16 
Z1=E+SIN  (2)  *  (2T*1  .0) 

?1=F+COS(Z)  *  (ST*1.6) 

GO  TO  103 

101  E1=i-*-SIl»  (Z)  *  (HT*2.2) 

E1=E+CO S  (2)  *  (8T*2.2) 

GO  TO  103 

102  E1=£+SIN  (2)  *  (KT*2.2) 

E1=E+COS  (2)  *  (5T*2.2) 

103  CONTINUE 

II  (Z-5. 235987)  202,205,205 

202  IF  (2-2. 7y2S2o4)  204,203,203 

203  XX=1.G 

GO  TO  25G 

204  II  (2-0.5235967) 205,221,221 

•  205  0=2.0 

GO  TO  250 

221  IF  (N-100.)  231  ,2ji4,234 

231  IE  (i.— 10.)  232,233,233 

232  0  =2.0 

GO  TO  250 

233  XX=3 . 0 

GO  TO  250 

234  XX=4.0 
250  CONTINUE 

TX=X+SIN  (0)  *  (5T*0) 

.  TY=Y+COS  (U)  MHT»XX) 
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C***»**»*  PLOT  SYEBQL 

122 

CALL  PLOT ( E ,F , 3) 

123 

CALL  PLOT  (X,I  ,2) 

124 

CALL  PLOT  (A, b, 3) 

125 

CALL  PLOT  (C, D ,2) 

126 

CALL  NUMBER  (El  ,F1  ,BT,DL,0  .0,-1) 

127 

I?  (CCNU3  -E0.  0.0) GO  TO  30 

128 

CALL  HUftBER  (TX,TY,HT,N,&.0,-1) 

129 

GO  TO  30 

130 

50 

CALL  SYJ1BOL  (X,T,2*B1 , 003,0. 0,-1) 

131 

30 

RETURN 

132 

END 

133 

C ********  SUBROUTINE  EESO  *»»* 

134 

C ********  TO  PLOT  TREND  AND  PLUttGE  Of 

135 

C********  CESOSCOPIC  FOLD  AILS 

1.36 

SUBROUTINE  fiESO  (X,Y ,1'R  ,PL  ,BT  ,IDES 

137 

REAL  X,Y,TR,PL,3T 

138 

INTEGER  IDENT 

139 

TE=ABS  (TR-450.) 

140 

TT=TB*1 . 7  4532  9E-2 

14  1 

B A=TT*  (3Q*1.745329E-2) 

142 

CA=TT-(30*1 .745329E-2) 

143 

If  (IDEBT-2)  600,600,700 

144 

C ********  CALCULATE  SLICKE3SIDE  SYMBOL 

145 

600 

E=X+5IN  (TT)  *  (3.1*3  .) 

146 

?=T  +  COS  (TT)  *(HT*3.) 

147 

E1=I  +  SIN  (TT)  *  (3T*5.) 

148 

?1=I*COS  (TT)  *  (HT*b.) 

149 

E2=X+SIB  (TT)  *  (3T*t>  .5) 

150 

F2=Y  +  COS  (TT)  *  (KT*6.5) 

151 

BAX--E1-SIR  (31)  *31 

152 

BAY=F1-COS (BA) *HT 

153 

CAX=E1-SIa  (CA)  *31 

154 

CAY=F 1— COS  (CA) »HT 

155 

GO  TO  a  00 

156 

C**«****»  CALCULATE  FOLD  AXIS  SI3B0L 

157 

700 

E-X+SIN  (TT)  *  (3T*3.) 

158 

P=Y+COS  (TT)  *(31*3.) 

159 

21=X+SIN  (TT)  *  (HT*tt.) 

160 

F1-T+COS  (TT)  *  (ST*8.) 

161 

Z2=X+S1R  (TT)  *  (31*  9  .5) 

1o2 

F2=I+COS(TT)  *  (HT*9 . 5) 

1©3 

BA  X—E 1-SIN  (EA)  *  (HT*1  .5) 

164 

3AI=F1— COS  (BA)  *  (HI*!. 5) 

1o5 

CAX-El-SIS  (CA)  *  (BT*1 .5) 

loo 

CAY-F  1-COS  (Ca)  *  (ST*  1-5) 

167 

£**♦♦**-**  PLOT  SYEBOLS  **** 

1o6 

800 

CALL  Pa. OT  (E,F,3) 

1«>  9 

CALL  PLOT  (El, FI ,2) 

170 

CALL  PLOT  (BAX, SAY  ,2) 

171 

CALL  PLOT  (CAX,CAI,3) 

172 

CALL  PLOT  (E1,?1,2) 

173 

CALL.NUKBZE  (I2,F2,iTT,PL,Ci.Q,-1) 

174 

RETURN 

175 

END 
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APPENDIX  6. 


Profile  plots  of  the  down-plunge  projections  of  the 
twenty-two  structural  domains  delineated  for  the  study  area 
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APPENDIX  7. 


Listings 
outputs,  used 
structure  cont 


of  the  programs,  with  appropriate  inputs  and 
to  generate  the  coordinates  of  points  for 
ouring . 


DATAFORM 


inputs  4=-PROFILE  (output  from  SECT5) 
5=DOMX 


outputs  6=PROFj.LE 


1  -  c 

2  C 

3  C 
u 

5 

5 

7 

S 

9 

10 

n 

12 

13 

14 

15 
1b 

17 

18 

19 

20  . 

F.tfD  Of  FI  Li. 


- L/.TA?Gf«K  IS  JjjS  IGMED  To  gifOSiV.T  PSGPILE  OUTPUT 

- —  FBOH  SECT  5  w  ITU  LATA  fSOE  DuRAIi;  FILES  IiiTO  A 
■ —  FORMAT  SUIT 6  nl,t  fOri  STRUCTURE  CONTOUalKG  TSOGHAnS 
IS  TIGER  D 
OO  1  1=  1  r  1  uO 

R  READ  (4, 10,m>-2)  I,  I,  FIT 
7  KC  A  D ( 5 , 1 1 ,  FM>-S )  o  ,  N  ,  A 
SD=D*1  . 

4  a’.S  ITS  (b  ,  1 2)  2  ,  l ,  f  IT,  SL>  ,  N  ,  .1 
GO  TO  6 

5  CC-KTIHUE 

H  B ITE  (6 , 1 3)  X  ,  \  ,  FI T 
i  COMIKUE 

1  CoSTIKOE 

.10  FORMAT  (bX,  3F12.  1) 

11  FORMAT (62X ,1- ,  0  X ,14 ,  2  X  ,  1 2  ) 

12  i’OKKAT  (4F8.1, 215) 

13  FORMAT (3Fb.1) 

2  Slop 
f'JD 


' 


- 
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HORIZON 


inputs  4=*SOOHC2* 

ddap  =  dip  direction  of  the  axial  plane 
diap  =  dip  of  the  axial  plane 
tnp  =  trend  of  the  normal  to  the  profile 
pnp  =  plunge  of  the  normal  to  the  profile 
cs  =  stratigraphic  position  (in  feet)  of  the 
coal  seam  above  the  formation  base 
xprof,  yprof  =  coordinates  of  the  lower  left- 
hand  corner  of  the  profile 
scale  =  scale  in  feet  per  inch 
5=PROFILE 

outputs  6=COALPROF  (file  containing  the  profile 
coordinates,  slope  and  weight  of  points 
calculated  to  be  on  the  coal  seam) 

9=profile  plot  of  the  points  on  the  coal  seam 


1 

C 

— 

hOallOH  IS  ofSIGMi*  TO  CAiCuLAiT  COOSiilii > TE5  A>JD 

2 

c 

— 

SLOP'S  ON  COAL  SKLl.  TO  ri  COMOUXrG  f.-iOT  PnOFILL 

3 

W 

— 

OUiCKOP  USIIxONS 

4 

1FTL32L  CS 

DT5 *1.74 5229 L-L 

6 

c 

7 

KILL  (4,  10)  LLA  i  ,Jlk£,Z::7  , ?.H«  ,  CS  .APLGf  ,  Xi?K05, SCALL 

S 

call  plots 

9 

c 

— 

DS.\W  50U1*’DaS1£S 

TO 

XP 9*1250^/10. 

11 

j-?a*-ypto?/ic. 

12 

SCA=SCALS/10. 

13 

CALL  11233  {1 .3,1  ,-1) 

1« 

CALL  AXIS2  (1.0, 1.0, ‘X-AXIS  *'iu*  , -11 ,  13  .0,0 .0, I?H,SCA  ,  1 

1b 

CALL  AXIS!  (11.0,1.0,*  •  , -  1 , 1 0 . 0 , 9 0 . 0 , VP? ,5CA 1 . 0) 

16 

CALL  IZIS 1 (1. 0, 1.0,  • t -L* IS  *  1  «.-■  * , 11, 1C. O, 90.0, y? 7  ,SCA  ,  - 

17 

CALL  AX1S2  (1.0,  11  .3  ,  »  *  ,  1 ,  1 0  .  0 , 0 . 0 , 1  IT. ,  5CA. ,  1  .  G ) 

Is 

TAF-D0LP**-1c30. 

19 

.-DIM 

20  . 

"i  ii  P  ^T  A?  — 4  sc  . 

21 

TSF=T!i?-450. 

22 

o::?-(T.*<P4  1s>  0. ) 

23 

ripa  (90  .—PNP)  *DT5 

24 

c 

— 

i»LTZXfil«E  DILSCTIOK  COS  IN  Is  OF  TL'3  XO^SAL 

25 

- — 

TO  iSL  AXIAL  PLAN L  And  IJ:  FOLD-AXIS 

26 

i  =  COS  (DTBn  ‘•.i  )  »CUS  (07 

27 

2«SIS  (DTK*lAi ) »CG£ (JTh»fAr) 

26 

CMlX  (1>I2*I-AI) 

29 

yaCC’S  (DTii* I  :;i  )  *CCS  (oiT.*PM) 

30 

--Sit;  (DTB*1  M  )  *  COS  {'J‘l 

31 

f*SUi  (t‘T2*rSi-j 

32 

— 

fi:-d  ti:s  Diu-cnos  cos  I  ms  of  i  '-rr  link  of  isi-sssrcuo*? 

33 

c 

) 

OF  T  22  ALi-iL  i-L/.M  AND  9  30  FI  Lb  (PiOJrCTIOM  LINS) 

34 

J  J 

SLa(’tf*F-C*4)/  (A*S-b*L) 

35 

6.1*  (-r-D*fi*,)  /£ 

3a 

Xf.=SCST  (1/(V:;L*»2+sr.**2)  ) 

37 

IL-k:i*bL 

3  6 

Y>  r.  =*X  t«  *  J-  m 

' 


. 

' 

■ 
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39 

c 

— - 

FISD  TSL  TLii-'D  iSL  Li.iU.oi  ,^r«  221  PITCH  0?  TUI 

40 

W 

— 

Ll:<£  OF  pioJlCHOfi  u *  li? 2  >r. OFILZ 

14*1  • 

2'J  1 

?IL=A  SI  H  (XI.) 

442 

TIL*DXF- (9v*aTi)  -ATAF  (COS  (DIP)  *TAH  (ASIF  (Si;,  (PILJ/SI.*;  (51?) 

<43 

D  IF* TIL  -  (D K l* -  9 0  *  D  IF. ) 

4<* 

A5P«-<iSI.i  (COS  (vlr)  *CCS  (50*013-? IL)  )  ) 

45 

AI*-4.o1/30G. 

14  n 

0 

w 

— 

A I  AD  IK  PSuFlLi  DATA 

47 

202 

uO  1  1=1,100 

43 

SiA5(5,  11,2::i>=2)  I,Y , PIT, SD 

49 

* 

rcT=?rr 

50 

S  LOPZ=- TA.4  (PILOTS) 

51 

IF  (PIT  .20.  -t.95.9)  GO  TO  100 

52 

IF  (PIT  .GT.  9 j . >  ? IT=? IT- 1 aO . 

53 

21u 

I?  (SO  .ZQ •  -995.5) GO  10  4 

54 

— 

ulCZXfilXZ  STiATlGAAPillC  S2?A5aTI0K  0 GTC.-;0? 

55 

w 

— 

a:;d  coal  szax 

55 

C0=CS-5D 

57 

CS=A££ (CD) 

53  ' 

c 

— 

calculate  V?-I5i:T 

55 

*r*2.7lb2t 15** (i.IvCh) 

to 

ruL*  (PIT-50.) *D7h 

61 

W 

— 

FI; to  ASGLL  sLS'tili  TS  aCZS  OF  AaIAL  FLiti*  5  aTL  aZL'DTKG 

62 

‘asg-app-pol 

63 

IF  (ANG  .GT.  7i*UTa)  GO  TO  5 

64 

1?\ASG  .  LT  .  •  -75*DTi)  GO  TO  5 

t>5 

CL=CD/COS  (V.'G) 

so 

** 

s* 

— 

FIND  PHOJLCTiO  CGOSoS  AK i)  SLOPS  OF  GGJLL  SSAf. 

67 

203 

X7  =  (X+COS  (AFP)  *CL) 

63 

i  1  -  (1  -Slii  {  AP?)  *CLJ 

69 

GO  TO  b 

70- 

5 

I7=2-*-COS  (PIT  1  *D T S)  *CD 

71 

Y1*r-*-SIH  (PII7*0TH)  *CD 

72 

e 

73 

o 

SPITZ  (o,12)  X7,il,SLOP2,  ST 

7  4 

GO  TO  5 

75 

700 

IF (ASS (CD)  .GT.  20.) GO  TO  4 

76 

SL0P8=0 .0 

77 

WT=1.U 

7a 

i»F ITZ  (6,1 2)  ;c,?,-lofl,wt 

79 

X  1=X 

30 

m*r 

31 

CO  ATI. '<0  2 

62 

X0=  (  (X7 SCiLi)  -ZP20F)  /SC A  LL 

S3 

10  =  (  (Tl^SCALL) -YPfOf) /SCALI 

gu 

. 

CALL  S1XS0L  (Xu,  10,0  .(,5,004,0.0,-7) 

85 

CONTISUI 

66 

7 

COG  Ti:<  n  z 

57 

10 

FOS^/T  (4FJ  .  7, 14 , 356 . 1) 

8b 

17 

F0l.fi  A'r  (hF5 . 1) 

69 

72 

FGr.ZA i  ( fb  .  1  , Fa  .  i  ,7 u  .4  ,P  10  •  o) 

90 

2 

Call  plot  ( c. .  r o . ,  9 59 ) 

91 

STOP 

92 

r  *n 

OF 

/Hi 

• 

' 
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SLOPE 


inputs  1=COALPBOF 
2=*SOUBCE* 

a  =  trend  of  fold-axis 
b  =  plunge  of  fold-axis 

xprof,  yprof  =  coordinates  of  the  lower  left- 
hand  corner  of  the  profile 

scale  =  scale  of  the  profile  in  feet  per  inch 
ri  =  distance  to  1%  influence  for  the 
weighting  function 

xmin  =  lowest  profile  X-coordinate  to  be 
calculated 

xmax  =  highest  profile  X-coordinate  to  be 
calculated 

3=DOMX  {with  additional  map  coordinates  added 
from  within  the  domain) 

outputs  5=*?BINT*  (listing  of  map  coordinates  and  slopes 
of  the  coal  seam) 

6=C0ALMAP  (map  coordinates  and  slopes  of  the  coal 
seam  for  contouring) 

9=COALPLOT  (profile  plot  of  the  coal  seam 
locations  projected  parallel  to  the  fold-axis) 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1C 

11 

12 

13 

14 

15 

1b 

17 

IB. 

19 

20 

21 

22 

2  J 

24 

20 

26 

27 

28 

29 

30 

31 

32 

•>.* 

ju 

30 

36 

37 
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C -  TO  »,LHi:i(AT>:  CUOLDI  !<  ATKS  for  structwkk  CONTOURING 

_  Dirt  LM SION  Thl  (9)  ,Ttt2  (5)  ,  XX  (200)  ,  I  i  (2U0)  ,  UP  (200)  ,  4T  (200) 

RKAL  Lr  ,rtP,  HI*  PM  ,  Jirt 

ISTKGhR  A  t'i ,  ?jI 
RKAL  OP. 

!( f:  A  LM  2 , 2  0 )  A  f  b  ,  X  I  h  <  •  P ,  V  I*  iUi  f ,  SC  A  L  £ ,  ft  I ,  X  W  N  ,  X  ft  A  X 

V H ITT,  (5,24)  i. ,  b 

A I -“4  .b 1/R J 

CALL  PLOTS 

XP!i*APbOF/lD. 

ir-u-YPi:of-/io. 
sc  a =sc  A  l  e/  i  o . 

CALL  AX2KP  (1.0 , 1  ,-1) 

CALL  AXJS2  (1.0, 1.0,  »X -AXIS  ♦ 1 0 • , - 1 1 , 1 0 . 0 , 0 . 0 , X PM , SC A , 1 .0) 
CALL  AXIS2  (11-0,  1.0, 1  1  ,  - 1 ,  10.0,90.0, 1‘  I*  ft  ,  SC  A  ,  “  1  .  0 ) 

CALL  A  X I S  2  (1 .0,1  .0,  *  Y-/.IIS  *10*,  11,  10.0,90.0,  YFR,SCA , -1 .0) 
CALL  AXIS 2  (1.0, 11.0, •  • , 1 ,  10 . 0 , 0 .0 , XWt , SCA , 1 .0) 

DT1<=*1  .7  4b  32 ‘»f -2 
A  — A  *D!fl 
B=0*£»Tft 
K- 1 

DO  1  J- 1,200 

HtfAD  (1,21,  hNI»=2)  XX  («J)  ,  i  i  ( J)  ,  J/P  (J)  ,W’L’  (J) 

M  =  R  ♦  1 

1  -COUT1  Mil  f. 

2  SA  =  SIN(a) 

CA  -COS  (A) 

SU  =  SlN(b) 

cn=cns(ii) 

C - CULATL  TK  ANhFlJhf'ATXu.i  F  ATh  LX  Tftl 

Tfl  1  (1)  =CA 
TP  1  (2)  -  -  S  A 
T r.  1  (3)  -0 . 0 
*J  U  t't)  aLA*LP 

■ml  {:*/ 

Th!  1(b)  -CL) 

Tfl  1  (7)  --S  A  *CU 


38 

39 

ao 

4  1 

42 

43 

44 

4  5 

4b 

47 

4b 

as 

50 

51 

52 

53 

54 

55 

56 

57 

53 

59 

60 

61 

62 

63 

64 

65 

66 

©7 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

31 

62 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 
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Tf*  1  (8)  =-CA*C  3 
.  T K  1  (9) =Sli 

C -  CREATE  TRANSFORMATION  RATS  17.  T»2 

7fc2(1)*-SA 
Id?  (2)  =CA*S 3 
Tf  2  (2)  =CA*C2 
TF.2  (4)  =CA 
7-2  (5)  -SA*S3 
7rt2  (6)  =SA*C2 
rr.2  (7)  =o.o 
T«2  (©)  -~C3 
M2  (5)  *SS 
►  SITS  (5,25) 

W  it  ITS  (5,26) 

1000  READ  0,22,  £XD=4)  X,  7,21 
TOP=Ll*  1.0 

C -  ERAKSFOEH  HAP  CG08DIN ATE5  TO  PROFILE  COORDINATES 

X?=X*TKl  (1)  (2) 

y  ?i=x*tki  (4)  +  Y*?r  i  (5) 

Z?=-  (X*Tfl1  (7)  +  Y*Tf5l  (8)  ) 

C - CALCULATE  TEE  PROFILE  COORDINATES  OF  7 ME  dOnIZOM 

1?  (XP  .LT.  XtlNJGu  TO  9 
IF  (XT  .GT.  XS  AX)  GO  TO  S 
S‘Xs0 . 0 
SDP=C.Q 

sw=c.o 

SH 0=0.0 
DO  100  1=1, S 
DX=XP-XX(I) 

Y  1=YY  (I)  +DX*uP  (I) 

DH=ABS  (DX) 

C  -  DETERMINE  WEIGHT 

'•*=(2.71828  1&*»  (Ai-DB)  )  *W T  (I)  * 
pfi  S  »'  =S  W  +  5 

SY=SY+Y 1 *4 
SD?=SD?+D?  (I)  *W 
IF  (DP  (I)  . EQ.  0.0)  GO  TO  1U1 

SVD=S»D+w 
101  CONTINUE 
100  CONTINUE 

??=sr/ss 

IF  (SVD  .EQ .  0.0)  GO  ZO  179 

C - DE7EB3ISE  SLOPE 

SLOP£=SDP/SWD 
GO  TO  180 
179  S1CPE=0 .0 


C 

n 

C 


-  C0E2ECT- COORDINATES  TO  SLOPE  OF  PROFILE 

180  YS2?=YP-T?1 

ZCHA=YSEP*TAH  (5) 

2P=2?-2CBA 

TO= ( (XP+SCALE) -XPriOF)  /SCALE 
YG  = ( (YpfSCALE) -YPROF)  /SCALE 
CALL  STHBOL(XO,Y 0,0 .05,004,0.0,-1) 

— TRANSFORM  PROFILE  COORDINATES  TO  EL?  COORDINATES 
XJ=XP*7E2  (1)  *Y?*Tr:2  (2)  +  ZP«TE2  (3) 

Y  J=Z?*T32  (4)  •*-Y?*Tfi2  (5)  -*»Z?*TN2  (6) 

ZJ  =-  (XP-T32  (7)  +  yp*T=2  (3)  ♦  ZP*T3  2  (9)  ) 

- DETE2EI2E  THE  APPARENT  DIP  OF  TEE  HOnliON  ALO.'LG 

-  HAP  AXES 

- ANGLE  BETWEEN  PROFILE  I  Y  AND  2  All 0  TRL  POLE  TO 


Of 


HI 

HE 


I NT  £  f  LIT 


son 


10  S 


- 


109 


56 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 
127 
129 
129 
120 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

1,4 

145 

146 

147 

148 
14  9  ' 

150 

151 

152 


C 


C 


2'J 


51 


201 


If  (SLOP E  .IT.  0 . 0)  GO  TO  50 
IP  (SLOPE  .  iy.  0.0)  GO  TO  500 
XA=ATAH  (1  ./SL0F2) 

Y  A  *—  (ATA:.  (SLOP  7.)  ) 

GO  TO  51 

XA=ATAH  (1 ./SLOPE) 

YA-ATAN  (SLOPS) 

COHTINOE 

UIHSCTI05  C0SIU2S  0?  THE  POLE 
LF*COS (90»DTk^YA) 

HP— CO  S  (S0*LTE+XA) 

N P=COS  (90* DTK) 

TltAHSFGEM  THE  DIEECIICN  COSIrii'S  TO  TS2  ML?  COORDINATE  SYSTEM 


LH=LP*TH2  (1)  -*BP*?B2  (2) 

EE* LP*TE2  (4)  +  rp*TE2  (5) 
r;3=—  (LP*TE2  {? ).  ♦BP  *TS2  (6)  ) 

CONVERT  THE  DIRECTION  COSXKIS  TO  A  .TBFKD  ARD  FLO  HOE 
IF  (HE  .LT.  -1.0)  NB=-1.0 

IF  (S3  .GT.  1.0)  PM=1.0 

PLUS— ASIK  (Kf.) 

TEN*  (LB/COS  (PLUi?)  ) 

IFfTBK  .LT.  -1.C)  TIUf*-1.0 
IF  (TEN  .GT  .  1.0)  Tif  N=1 .0 
TRSS=AFIR  (Th>) 

CGK1SBT  TLE  7EFNP  AMD  PL3RGE  TO  k  DIP  DIFSCTIOsi  ASD  DIP 
DX rDI5= ( 1 80*  DTK)  +TSEH 
Dir*  (90*DT1)  -PLKl* 

FIND  TSE  APPAHTNT  DIP  OF  il S?  EOSIZOii  OK  TEE  R-S  (X)  AXIS 

OK  THE  i-W  ( Y)  2.IIS 
ADX=SIR  (DIIDI K)  *TAX  (LLP) 

ADY-SIN  (S0*I»TE— DIPUIH)  *T  A  N  (DIP) 

GO  TO  501 


A  KD 


-  CALCULATE  OfEHb  HEXES*  THICK K  i-SS 

Sue  ADX  =0 . 0 
ADY=0.0 
5C1  GS  — T0P-2J 


-  WHITE  NE4  CUGRD-IR  ATF5 

WHITE (5,23) Y J , X J , 2J , ADA, ADX, TO?, 05 
WHITE  (6  ,23)  YJ  ,XJ  ,  SJ  ,  ADX  ,  ADY  ,  I  OP,  03 
9  CON'TI HUE 
3  CONTINUE 
GO  TO  1000 

20  FORMAT (8F6 . 1) 

21  FOSHAT  (2?6 . 1 ,  Po  .  4  ,F10 . 6) 

22  FORMAT (IX, 315) 

23  FORMAT  (7F 1 1  .3) 

24  FOEEAT ( *  1  *  ,5X , 1 GHFOLD  AXIS  ,2rF.1//) 

25  FOEEAT  (5X,21EP01.KTS  FOR  CON’TOU  2IKG/) 

2o  FOEEAT  (7X  ,  •  X » ,  1  u  X,  *  Y  •  ,  10X  ,  •  Z  »  ,  7X  ,  »X -SLOPE  •, 
*4X, ' I -SLOPE  *  ,  &X  ,  'TOPG’,77, *05') 

27  FORMAT  (  ‘ENTER  THE  IBTEE7AL  OF  IT  I-IFLCL’MCE’ ) 
2(>  FORMAT  (F8 . 1) 
u  CALL  PLOT  (C  .  ,0.  ,999) 

STOP 
SR  D 


. 


. 


APPENDIX  8. 

The  following  composite  stratigraphic  section  of  the 
Lascar  Formation  was  compiled  from  partial  sections  measured 
on  Packrat  Creek,  between  grid  points  28700E  23600N  and 
27230E  24480N,  Campbell  Flats,  from  27230E  24480N  to  26300E 
2677 ON,  Hells  Creek,  from  25350E  25950N  to  23200E  25900N, 
and  on  Grande  Mountain  from  37750E  11050N  to  37850E  10200N. 


Feet 

M 

Unit  D 

1898.5 

578.6 

Interbedded  sandstone,  shale  and  siltstone; 
abrupt  contact  with  overlying  Shaftesbury 

Formation 

1472.5 

448.8 

Granular  and  pebbly  sandstone 

1452.5 

442.7 

Medium-grained  sandstone 

1412.5 

430.5 

Granular  and  pebbly  sandstone 

1387.5 

422.9 

Interbedded  sandstone,  mudstone  and  siltstone 
with  a  4  ft  (1.2  m)  pebble  conglomerate  at  top 

Unit  C 

1270-5  387.2 

Eidge-f orming,  lithic,  medium-  to  fine-  grained 
sandstone;  flaggy  in  places  with  some  siltstone 
and  mudstone  bands 

1200.5  365.9 


1120.5 

Interbedded  mudstone  ,  siltstone  and  sandstone 

341.5 

Carbonaceous  mudstone 

1110.5  338.5 

Sandstone 


•  . 
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1107.5 

337.5 

Mudstone 

1104.5 

336.6 

Slightly  carbonaceous,  cross-bedded,  medium- 
grained  lithic  sandstone 

1095.5 

333.9 

Mudstone 

1081.5 

329.6 

Coal,  #11  seam 

1077.5 

328.4 

Carbonaceous  shale 

1076.5 

328.1 

Silty  mudstone  with  some  siltstone  bands 

1062.5 

323.8 

Medium-grained,  lithic,  poorly  indurated,  ridge¬ 
forming  sandstone  with  large,  oval,  reddish- 
brown,  limonite  stains  on  weathered  surfaces 

980.5 

298.8 

Shale  with  clay  ironstone  nodules,  #10  coal  seam 
is  covered  but  should  lie  between  950  and  980  ft 
(289. 5  and  298.  7  m) 

965.5 

294.3 

interbedded  fine-  to  medium-grained  recessive 
sandstone,  siltstone  and  some  carbonaceous  shale 

949.5 

289.4 

Very  fine-grained  sandstone 

947.5 

288.8 

Very  f ossilif erous,  orange-brown  weathering 
siltstone 

943.5 

287.6 

Mudstone  with  siltstone  bands 

939.5 

286.3 

Siltstone  with  some  sandy  bands 

925.5 

917.5 

282.1 

Shale 

279.6 

flaggy,  medium-grained,  cross-bedded,  slightly 
carbonaceous  sandstone  with  some  reddish-brown 
weathering  clay  ironstone  nodules 

917,5 
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899.5  274.2 

Shale  with  some  sandy  and  silty  bands 

876.5  267.1 

Coal  and  coaly  shale 

871.5  265.6 

Medium-grained  argillaceous  sandstone 

869.5  265 

Interbedded  shale  and  siltstone 

840.5  256.2 

#8  seam 

835.5  254.6 

Shale  with  some  fine-  to  medium-grained  sandy 
zones  at  top  and  bottom 

812.  5  247.6 

Carconaceous  and  coaly  shale 

807.5  246.1 

Siltstone  with  abundant  plant  fossils 

802.  5  244.6 

Carbonaceous  shale 

791.5  241.2 

Coal 

786.5  239.7 

Fine-grained  silty  sandstone  with  abundant  plant 
fossils 

761.5  232.1 

Mudstone 

756.5  230.6 

Fine-grained  silty  sandstone 

741.5  226 

Mudstone  with  some  coal  from  727  to  729  ft  (221.6 
to  222.2  m) 

709.5  216.2 

Coal  and  coaly  shale 

706.5  215.3 

Light  grey  weathering,  blocky  siltstone 

698.5  212.9 

Bidge-f orming ,  massive,  fine-grained,  silty 


. 


sandstone,  locally  referred  to  as  the  ’Super  4 
Sandstone • 


638.5 

629.5 

618 

607 

599 

585 

550 

545 

487 

459 

458 

421 

401 


194.6 

Mudstone,  massive  near  top,  bedded  near  base 

191.9 

Coal,  #4  seam 
188.4 

Fine-grained,  argillaceous  sandstone;  silty  in 
the  upper  2  ft  (0.6  m) 

185 

Mudstone  with  silty  bands 

182.6 

Fine-grained  sandstone 

178.3 

Mudstone,  silty  in  the  upper  10  ft  (3m); 
contains  the  ’Clam  Zone* 

167.6 

Coal,  #3  seam 
166.1 

Medium-grained  lithic  sandstone  which  weathers  a 
light  grey  with  some  reddish-brow n  limonite 
staining 

148.  4 

Interbedded  mudstone,  laminated  siltstone  and 
thinly  bedded  fine-grained  sandstone 

139.9 

Chert  pebble  conglomerate  with  medium-grained 
sandstone  matrix 

139.6 

Mudstone  and  siltstone  with  sandstone  bands  0-5 
to  1.5  ft  (0.2  to  0.5  m)  thick 

128.3 

Fine-grained  sandstone  with  silty  and  shaly  band 


Unit  B 

122.2 

Mudstone  with  rare,  orange-brown  weathering, 
nodular  ironstone  bands  and  rare  thin  bedded  and 


. 


' 


cross-bedded  fine-grained  sandstone  bands 


Unit  A 

299.5  51.3 

Carbonaceous,  argillaceous  and  sandy  laminated 
siltstone  with  abundant  chert  pebbles 

295.5  90.1 

Partially  covered  fine-grained  sandstone, 
carbonaceous  near  the  base 

275.  5  84 

Carnonaceous  siltstone  and  shale 

271.5  82.7 

Fine-grained  argillaceous  sandstone 

262.5  80 

Interbedded  shale,  siltstone  and  laminated 
siltstone  with  some  cross- bedded ,  fine-grained 
sandstone  bands 

245.  5  74.8 

Covered 

219.5  66.9 

Massive  fine-grained  sandstone  with  rare  bands  of 
medium-  and  coarse-grained  sandstone  and  some 
lenticular  and  cross-bedded  zones 

198.5  60.5 

Mudstone  becoming  silty  and  sandy  near  the  top 

190.5  58.1 

Thin  bedded,  fine-grained,  carbonaceous  sandstone 
which  grades  into  silty  sandstone  near  the  top 

176.5  53.8 

Slightly  sandy  siltstone  with  abundant  plant 
fossils 

168.5  51.4 

Coarse-  and  medium-grained  sandstone  with  some 
interbedded  fine-grained  sandstone  in  the  upper 
half,  pebbly  near  the  base 

148.5  45.3 

Chert  pebole  conglomerate  with  medium-  and  coarse 
grained  sandstone  matrix 


■ 


115 


146.  5 

144 

142.  5 

137.5 

132.5 

129.5 

128 

126 

124 

119 

114 

107 

101 

85 

70 

66 


44.7 

Silty  shale 

43.9 

Fine-grained  argillaceous  sandstone 

43.4 

Hard,  carbonaceous  mudstone  with  abundant  plant 
fossils 

41.9 

Thin  bedded,  fine-grained  sandstone  and  siltstone 

40.4 
Shale 

39.5 
Coal 

39 

Shale 

38.4 

Beddish-brown  weathering  siltstone 

37.8 

Covered,  appears  to  be  mainly  shale  but  there  may 
be  a  thin  coal  seam  present 

36.3 

Medium  grey  weathering,  fine-  to  medium-grained 
argillaceous  sandstone 

34.7 

Medium-grained,  brownish  weathering,  cross- 
bedded,  flaggy  sandstone 

32.6 

Partially  covered,  nodular,  reddish-brown 
weathering  mudstone 

30.8 

Coarse-and  medium-grained  sandstone 

25.9 

Coarse-  to  medium-grained,  pebbly  sandstone  with 
common  lenticular  pebble  conglomerate  beds 

21.3 

Coal,  #1  seam 

20.  1 

Interbedded  shale  and  fine-grained  sandstone  with 


<s*=% 


rare  medium-grained  sandstone  bands 


58.5  17.8 

Nodular,  silty,  reddish-brown  weathering  mudstone 

57.5  17.5 

Shale  with  some  thin  fine-grained  sandstone 
bands,  plant  remains  and  carbonaceous  and  coaly 
zones 

47.5  14.5 

Silty  mudstone  with  abundant  fossil  plants 

45.5  13.9 

fine-grained,  silty  sandstone 

44.5  13.6 

Carbonaceous  shale  with  numerous  plant  fossils 

31.5  9.6 

Fine-grained,  silty  sandstone 
29  8.8 

Carbonaceous  shale  with  abundant  plant  remains 
27  8.2 

Reddish-brown  weathering  siltstone  with  some  clay 
ironstone  nodules 

26  7.9 

Interbedded  siltstone  and  silty  sandstone  with 
some  fossil  plant  remains 

0 

Cadomin  Formation 
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